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1. INTRODUCTION 


The Applied Physics Laboratory (APL) has conducted 
an experiment for the U.S. Air Force (USAF) and _ the 
Defense Advanced Research Projects Agency (DARPA) to 
determine the validity of one facet of the theory of 
groundwave propagation at 100 kHz. The theory indicates 
that, if proper account is taken of group and phase 
velocities in the propagation of pulses of the Loran-C 
radionavigation service, geodetic position should be 
computable using the vacuum speed of light and_ the 
times of arrival of the pulse and the carrier. 
Specifically, the test goal was to determine if an 
analytic function can be developed for operational use 
that relates secondary phase factor (SPF) to 
envelope-to-cycle difference (ECD) so that geodetic 
position can be computed accurately and in real time. 


Test data were obtained in the eastern United 
States using the East Coast Loran-C chain and _ two 
identical measurement systems, one in a mobile test van 
and the other at a fixed site at the U.S. Naval 
Observatory (NAVOBSY). The van occupied 10 field sites 
during the course of the test. APL was assisted in 
this test by NAVOBSY and the Defense Mapping Agency 
Topographic Center (DMATC). The positions of all sites 
were precisely surveyed by DMATC using the Navy 
Navigation Satellite System. The sites are on ray 
paths from Loran-C stations through NAVOBSY. A’ team 
from NAVOBSY executed time transfers to each field site 
during the collection of loran data. The collected 
data include very accurate measurements of the Loran-C 
signal voltage and the absolute time at which they were 
made at each of the 10 pairs of sites. 


The test operations, including site selection and 
surveys, acquisition of Loran-C and time transfer data, 
and a Loran-C receiver response test, were initiated in 
April 1974 and completed in November 1975. This report 
(Volume D of TG 1298) documents the data, its analysis, 
and findings. Volume A is the summary report, Volume B 
describes the test operations, and Volume C describes 


the measurement system. All volumes are written under 
the assumption that the reader is reasonably familiar 
with the Loran-C system of radionavigation. Reference 


1 is suggested reading for those who are not. 


As a result of extensive data analysis, it was 
possible to establish statistically significant 
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measures £ SPF and ECD for each of the 12 tests 
conducted at the field sites and at NAVOBSY. The 
voltage measurements were used also to isolate the 
modulations of the carrier that give the pulses from 
the Carolina Beach, Dana, and Nantucket Stations their 
unique characteristics. The pulses from each station 
contain both in-phase and in-quadrature modulations of ; 
the carrier that make it impossible to make direct 
observations of the carrier. Largely because of this, 
the functional relationship between ECD and SPF that, 
it was hoped, would be derived from the data is 
presently not well defined. However, trends are 
observed in the data that tend to support the 
hypothesis that the desired functional relationship 
exists, at least under certain conditions. Additional 
analysis and perhaps’ tests will be necessary to define 
these conditions or to prove conclusively the existence 
of such a relationship. 
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2. OVERVIEW 


Volume D of TG 1298 describes the processes that 
were applied to the portion of the field data that 


pertains to the primary goal of the experiment. rt 
also records the processed data, a number of 
intermediate results, and findings regarding the 


attainment of the goal. We have chosen not to publish 
the raw field data and the rather voluminous computer 
output that was generated in developing an 
understanding of the data and their statistical 
characteristics. 


The field data are recorded on approximately 100 
10-in. reels of magnetic tape and will be stored 
indefinitely. Some of the reels of tape contain 
recordings of ancillary data. From one set of these 
data, the time difference coordinates of the satellite 
survey monuments at the various field sites could be 
determined. A second set was recorded to quantify any 

§ diurnal effect on Loran-C emissions. The ancillary 
data have not been processed, and at the present time 
there are no plans to do so. 


The complete report is in four volumes. Volume A 
is a summary report that touches on all aspects of the 
experiment. It includes excerpts from this volume and 
discusses the findings of the experiment. Volume B 
describes test operations, and Volume c is a 
description of the measurement systems used to gather 
data. Volume D_ starts at the point in the experiment 
where all the observations of the Loran-C’ signals, 

2 along with other pertinent information, are available 
4 as recordings on magnetic tape. 


If the reader wishes to understand this volume in 
relation to the complete experiment, it is suggested 
that Volume A should be read f£i¢se. Otherwise, the 
reader will be exposed only to the collected data, its 
processing, and analysis. 
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RECORDED DATA 


The data recorded on the magnetic tape are 
identified in Table 1. These data were recorded in 
blocks that were limited in size to 512 8-bit words 
(bytes). The format within each block is shown in Fig. 
1 


| The analog to digital (A to DP) converters were 
I adjusted such that 1 V produced the largest binary 
output number, i.e., 2 to the 13th power minus 1. Since 
each voltage measurement consisted of 13 bits plus a 
| sign bit, two bytes were required for each measurement. 


The locations of the voltage measurements on the 


Loran-C pulse are shown in Fig. 2. The odd-numbered 
Measurements are shown as "x" bits and the even as "y" 
bite In Hig. +. The block size accommodated all the 


voltage measurements on seven pulses, the Universal 
Time Coordinated (UTC) time of the first x sample in 
the block, and all the "header" data read from the 
Interface Unit (IU) of the measurement’ system. Four 
hundred seventy-seven bytes were used in each block, 


The locations of the field sites and the days on } 
which data were taken are given in Table 2. The exact 
geodetic positions of the survey monuments over which 
the antenna was placed at each field site and NAVOBSY 
are given in Volume B. AG. @8ch Field - site and 
simultaneously at NAVOBSY, primary data were recorded 
on the magnetic tape in three files, each containing a 
20-min segment of the data. The first file was started 
a6 330) PM; the second at Soh PM; ‘and the third et 
42b2 PM -EDE« No changes were made among recordings 
other than to indicate, if necessary, a change in the 
temperature of other environmental conditions on the IU 
control panel. Recordings of primary data were made 
with the voltage sample points located in time as shown 
in Fig. 2, with two exceptions. At Marietta, OH, the 
back side of the antenna was inadvertently pointed 
toward the Dana Station, and the voltage measurements 
were made 5 us later than shown in Fig. Ze Later 
measurements made specifically for the purpose (5 us 
earlier than shown in Fig. 2) showed that the effect of 
reverse pointing was limited to an inversion of the 
phase, 
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Table 1 


Taped Data 


Parameter Digits/Bits Units 

Voltage at Odd Samples, x 13 bits + sign V 
Voltage at Even Samples, y 13 ‘bits + sign V 
Attenuation 2 digits dB -7 
Time 13 digits hs min, 10 s 
Day Number 3 digits - 
Year 1 digit se 
Group Repetition Interval 3 digits 10 s 
Station Code | digit - 
Site Code 2 digits - 
Weather Code 2 digits we 
Atmospheric Pressure 5 digits 10) in. Hg 
Wet Bulb Temperature 3 digits 10_,°¢ 
Dry Bulb Temperature (outside) 3 digits 10_,°¢ 
Dry Bulb Temperature (inside) 3 digits Lo YE 
Data Recording Mode 1 digit ee 
Predicted Time Calibration 9 digits + sign L@ @ 
Note Code 1 digit = 

€ 
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Byte (most/least significant bit) 
1 
2 
3 (112 pairs of x, y samples, 7 GRI's, 
4 in 2's complement binary) 
448 
449 
Minute 
Second Second 
Second Second UTC time of first x sample in block) 
Second Second 
Second 3 Second 
455 Second F wer 
_ Atten. (10!) 
Atten. PEC (+) 
PTC (10-2) 
PTC (10-4) (PTC = predicted time calibration) 
460 - PTC (10-8) 
PTC (10°8) 
Day (102) 
Day (109) 
Xmitter 
Site (101) 
Wx (101) (Wx = weather code) 
Note 
DBT! (100) (DBTI = inside dry bulb temperature) 
WBT (101) (WBT = wet bulb temperature) 
470 WBT (10-1) (DBT = outside dry bulb temperature) 
DBT (109) 
BARO (101) 
BARO (10°!) (BARO = barometric pressure) 
BARO- (10-3) 
475 GRI (10-1) 
GRI (10°3) 
-- = zero-filled characters) 
512 


Fig. 1 Raw Data Tape Format for Each Block of Recording 
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Site 


Nantucket 1 


Baseline 
Dana 3 
Dana 2 
Dana 1 


Nantucket 3 
Nantucket 2 


Carolina 
Carolina 
Carolina 
Carolina 
Baseline 


wrreh 
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Table 2 


Location 


Toms River, NJ 
NAVOBSY, DC 
Georgetown, DE 
Marietta, OH 
Danville, IN 
Bluefield, WV 
Grottoes, VA 
Emporia, VA 
Wilmington, NC 
Dexter, NY 
Towanda, PA 
NAVOBSY, DC 
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Test Field Sites Visited 


Test Date (1975) 


29 
3 
5 

10 

12 

a 

19 

24 

26 

390 
2 
7 


May 

June 
June 
June 
June 
June 
June 
June 
June 
June 
July 
July 
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During each recording session at each field site 
and at NAVOBSY, calibration data were also recorded in 
separate files on the tape. These data showed that the 
measurement systems retained their calibrations 
throughout the test. The calibration method is 
discussed in Volume B, 


DATA PROCESSING 


Refer to Fig. 3. The format of the data recorded 
on the tapes labeled "Field Site" and “Fixed Site" is 
shown in Fig. 1. The data from both of these tapes, 
together with job control cards, were subjected to 
preprocessing immediately upon the return of the tapes 
to the Laboratory on data-taking days (see Volume BB). 
The preprocessor disregarded all voltage measurements 
over 1 V (first-stage editing). Its printed output 
consisted of the means and standard deviations of 
groups of 50 contiguous pulses and the grand means and 
standard deviations of the voltage measurements on all 
the pulses accepted from each 20-min segment of data. 
The data taken during calibration recordings were 
processed in the same way as_ those for pulses. The 
printed output of the preprocessor was used to assess 
the quality of the data prior to moving to the next 
site. The output also included a condensed data tape 
containing these grand means and standard deviations, 
as well as all the acceptable pulse data in a revised 
format. A deck of punched cards was also produced that 
recorded an editing criterion, taken to be twice the 
average of the standard deviations at each sample point 
along the pulse for use in further editing of the data. 
The second-stage editor rejected all the data for any 
pulse if any voltage sample of that pulse deviated from 
the mean by more than the editing criterion. The 
printed output of the second-stage editor was the same 
as that out of the preprocessor. The punched card 
output recorded the means of the accepted voltage 
samples taken during each 20-min data segment. 


The processing shown in Fig. 3 up to the tape of 
condensed data was done in two stages on the IBM 360/91 


computer. In the first stage, the first 16 voltage 
samples of the pulse were processed; in the’ second 
stage, the last 16. The computer programs for the data 


processing are documented in Appendix A. 

There is one reel of tape of condensed data for 
each site and one for each recording of baseline data. 
There are 12 files of data on each of these tapes. 
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Files 1 through 6 contain all the edited data for 
sample points 1 through 16, and files 7 through 12 for 
sample points 17 through 32. Files 1, 2, and 3 and 7, 
8, and 9 contain 20-min segments of data recorded at 
NAVOBSY: files 4, 5, and 6 and 10, Il, and 12 contain 
the simultaneous recordings at the field site. Each 
file is organized into a succession of records. The 
first record is the header data, i.e., the data 
recorded from the front panel of the IU (see Volumes B 
and C). There are as many more records as there are 
pulses that survived first-stage editing. Each of 
these records contains the following data in the order 
mentioned: the UTC of the first voltage sample in the 
record, 16 mean voltages, 16 standard deviations of 
these voltages, and 800 voltage measurements, i.e., 16 
Measurements on 50 pulses. The preprocessor, in 
organizing the condensed data, correlated the time of 
sampling each pulse so that samples with the same index 
number on corresponding records for NAVOBSY and the 
field site are samples of the same pulse. This 
requires that both records be of the same length. This 
was done by making the field site data record the same 
length as NAVOBSY, either by truncating it or filling 
it with obviously false data (flags). Subsequently, 
the second-stage editor recognized these flags and 
processed only those data for which there were _ two 
observations of the same pulse. 


ANALYSIS 
Reduced Data 


Analysis of the experimental results was performed 
on the mean values of voltage measured during each 
20-min segment of data. These are the mean voltages 
recorded on the punched card output of the second-stage 
editor (Fig. 3). Appendix B contains a table of these 


data (Table B-1l). Also included in Table B-l are 
"Other Data," namely, the date, the time of the first 
data record in each 20-min segment, the difference in 


UTC time between NAVOBSY and the field site, the 
distance from the transmitting station to the recording 
site, and the attenuation setting of the attenuator 


(see Volume B). Appendix B also includes a table of 
the signal-to-noise ratios computed on ae per-pulse 
basis for the first segment of data. These 


signal-to-noise ratios are indicative of the quality of 
edited data. 
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Polynomial Fits 


The first step in the analysis was polynomial 


fitting of the mean voltages. The times of the voltage 
Measurements along the pulse were very accurately timed 
at 2.5- and 7.5-us intervals in synchronism with the 
output of the cesium standards. The odd=- and 
even-numbered voltage measurements, therefore, are 
exactly in quadrature to each other when sampling a 
100-kHz signal. To sample a loran pulse, the 
odd-numbered samples were taken near zero crossings, 
and the even near the peaks, Ls@ss in phase. In 


anticipation of the advantages of such orthogonality, 
the mean in-phase (I) voltages were fitted separately 
from the mean quadrature (Q) voltages. Although most 
difficult to confirm experimentally, it is assumed that 
the loran pulse starts with zero voltage and zero 
slope. To accommodate this assumption, it was decided 
to fit the I and Q mean voltages with polynomials with 
the coefficients of the first two terms (ascending 
order of degree) set to zero. After considerable 
experimentation, it was found that 10th degree 
polynomials fitted the mean voltages such that the 
maximum residuals were the same order of magnitude as 
the standard deviation of the voltages themselves. 
Because of the constraint imposed on the first two 
coefficients, these fittings required an iterative 
procedure to determine the difference in time between 
the time reference of the voltage samples and the start 
of the polynomials. Then, to realign the sample points 
in time, the polynomials were translated by this 
difference in time. The first two coefficients of the 
resulting polynomials are no longer zero. These 
polynomials were judged to be adequate representations 
of the measured data for use in later illustrations of 


the pulses and subsequent analyses. The computer 
program used for the polynomial fits is listed in 
Appendix C. The coefficients of all the fitted 


polynomials are presented in Appendix D. 
Carrier Modulations and Pulse Shape 


Initial plots of the I and Q polynomials indicated 
that the quadrature component was quite substantial in 
all of the observations. Since this could have been 
the result of misalignment between the zero crossings 
in the pulse and the Q (odd-numbered) sample points, 
the I and Q polynomials were jointly manipulated to 
determine whether or not there was a position on the 
pulse for the sample points where all the Q voltages 
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would be zero. This was done by developing an analytic 
solution for the sample point locations for which the 
energy in the quadrature component of the pulse was a 
minimum. The algebraic development of this solution, 
as well as the computer program to execute it, are 
described in Appendix C. 


Appendix E contains plots of the 1 and Q 
components of the signal after adjustment for minimum 
energy in Q between 10 and 70 us. Since the I and Q 
components are orthogonal at the 100-kHz frequency of 
the Loran-C carrier, the I component can be considered 
to be the modulation of the sine of the carrier, and Q 
the modulation of the cosine of the carrier; that is, 


V(t-t0)=I(t-t0) sin(w(t-t0)+4)+Q(t-t0)cos(w(t-t0)+¢), (1) 


where V is the signal voltage, 
t is the time, 
td is the reference time, 


I(t-t0) is the I polynomial, 
Q(t-t0) is the Q polynomial, 


w is the angular frequency 
corresponding to 100 kHz, and 
(oy is the phase angle measured at t0. 


The shape of the pulse envelope was computed as q 
the square root of the sum of the squares of voltages § 
evaluated from the I and Q polynomials (see Appendix E 
C). The envelope shapes are shown in Appendix E. 


The Correction Function 


The envelope shapes of Appendix E are presented in 
such a way that a comparison can be made between the 
shapes observed at the field sites and the shapes 
observed simultaneously at NAVOBSY. The purpose of 
these comparisons was to disclose the differences 
between observations of the pulses at the field sites 
and at NAVOBSY, and to quantify the difference in pulse 
arrival times by one number, which we called ECM 
(envelope-to-cycle measure). 


An inspection of the envelope shapes presented in 
Appendix E indicates that ECM is many-valued; that is, 
it depends on the location on the pulse where the time 


comparison is made. This is a consequence of very 
complicated propagation effects acting on both the sine 
and cosine components of the signal. It was decided, 


therefore, that ECM would be evaluated as it would be 
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in a receiver that was perfectly phase tracking the 
Loran-C signal at the standard track point, that is, at 
the zero crossing at the end of the third cycle in the 
pulse, and also had the capability of measuring the 
envelope voltage at the same time as the phase track 
point. 


The inferred locations of the zero crossings 
computed using Eq. (1) are presented in Tables F-1 and 


F-2 in Appendix F, Note that none of the zero 
crossings are 5 wus apart, as they would be if the 
carrier were only amplitude modulated. Note also in 


these tables that the zero crossing at the end of the 
third cycle is not at 30 ws under the condition of 
minimum energy in Q(t) over the leading edge. To align 
the standard track point of all the observed pulses at 
30 ws, energy was shifted from I(t) to Q(t) such that 


the fractional parts of the "30.xxx-ys" columns of 
Tables F-l and F-2 were removed (to about three decimal 
places). Using these revised polynomials for I(t) and 


Q(t), the zero crossings were recomputed using Eq. (1), 
and the pulse envelope shapes were also recomputed (see 
Appendix C for computing programs). These zero 
crossings are presented in Tables F-3 and F-4. 
Appendix F also presents expanded views of the leading 
edges of the envelopes as they appear, having aligned 
the standard track points at 30 us. 


An inspection of Figs. F-1 through F-12 indicates 
that the repositioning of the zero crossings does not 
resolve the problem of multivalued ECM. We continued, 
therefore, to approach the problem as our hypothetical 
receiver would. ECM was computed relative to NAVOBSY 
by computing the voltage difference between the field 
site and NAVOBSY observations and converting this 
voltage to time using the average slope of the _ two 
envelopes at 30 us. The resulting values of ECM are 
presented in Table F-3. 


SPF has usually been defined as it pertains to the 
propagation On a 100-kHz continuous wave (CW) 
groundwave. Since the 100-kHz carrier cannot be 
observed in pulses when they are transmitted with 
cosine modulation, this definition lacks application to 
the problem at hand (see Volume A for further 
discussion on this point). Therefore, we are faced 
with anew concept, namely, a secondary phase measure 
(SPM) defined with respect to the arrival time of the 
standard third-cycle zero crossing as follows: 
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SPM = (D/C) = T3, (2) 


where D is the distance over which the loran 
pulse propagates, 
C is the vacuum speed of light, and 
T3 is the time of arrival of the zero 
crossing at the end of the third cycle. 


| 
| 
| 
The vacuum speed of light used in computations was 
299 792 456.2 m/s. The distance, D, and the time of 
arrival, T3, are computed from the basic data relative 
to their values that apply at NAVOBSY. The values of D 
and T3 used in Eq. (2) and the resulting values of SPM 
are included in Table F-3 of Appendix F. The value of 
T3 was computed as follows: 
| 


T3 = BIV = BIO = AT + AAT, (3) 


where BTV is the block time at the field 

site, i.e., the time of the first 
voltage sample in a 20-min data 
segment, 

BTO is the block time at NAVOBSY, 

AT is the time transfer result, and 

AAT is the difference in time adjustments, 
computed relative to NAVOBSY, 
required to align all the standard | 
third-cycle tracking points at 30 us. 


The values of BTV, BTO, and AT are tabulated as part of 
Table B-l of Appendix B. The value of AAT is obtained 
as an output of the computer programs used to align the 
standard third-cycle zero crossings. 


| 
The values computed for SPM are plotted against ) 

ECM in Fig. F-13 of Appendix F and in Fig. 26 of Volume | 

A. This figure is a representation of the correction 

function that was sought as the principal goal of the 

experiment. SPM is plotted against distance from 

NAVOBSY in Fig. F-14, and ECM against distance from 

NAVOBSY in Fig. F-15. The attenuation of the signal is i 

plotted against distance from the transmitters in Fig. 

F-16 of Appendix F and in Fig. 27 of Volume A. The 

significance of Figs. F-13 and F-16 is discussed in 

Volume A, 
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Frequency Filtering 


As noted in Volume A, Fig. F-13 tends to describe 
a low-order functional relationship between SPM and ECM 
for only four of the twelve possible triads of points 
plotted in Fig. F-13. These four are connected by a 
straight line in Fig, F=E3, In an attempt to 
understand why the remaining points appear to be in 
disarray, histograms were plotted initially of the data 
taken at NAVOBSY and simultaneously at Danville, IN, 


during the first 20-min recording. These histograms 
are shown in Appendix G as Figs. G-l and G-2, 
respectively. Contours of equal frequency of 


occurrence are reasonably circular and concentric in 
Fig. G-2, indicating that the dominant interference is 


randon,. The marginal distributions are reasonably 
normal. This is not true in Fig. G-l. The inclined 
barbeli contours indicate the presence of an 


interference that bears a coherent relationship with 
the loran signal. The distribution is double humped in 
spite of the fact that the marginal distributions are 
reasonably normal. 


The recorded data were then subjected to a Fourier 
analysis to identify the coherent interference, if any, 
and to subsequently remove it. The results of the 
Fourier analysis are shown in Figs. G-3 through G-6 in 
Appendix G. Each point plotted in these figures is the 
sum of the squares of 16 adjacent Fourier coefficients. 
The plots, therefore, display the power’ spectral 
density at reduced frequency resolution. Figure G-3 is 
the power spectral density of the quadrature data, and 
Fig. G-4, the in-phase data, taken at NAVOBSY. Figures 
G-5 and G-6 show the power spectral densities for the 
simultaneous Danville, EN, data. Note that at 
Danville, IN, there is an interference at frequency 
number 124 and at NAVOBSY at frequency number 154. 
These interferences were removed by detecting those 
frequencies observable in the plots that had power 
levels that exceeded K times the average power level of 


all frequencies. Effective filtering was accomplished 
with K=128. Then all 16 Fourter coefficients 
associated with each detected frequency and 32 


additional coefficients, i.e., 16 on each side, were 
set to zero. The resulting sets of coefficients were 
then subjected to the inverse Fourier transform, which 
reconstitutes the data set. Figures G-7 and G~-8 show 
the histograms of the frequency filtered data. 
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Because of the difficulties with the Georgetown, 

DE, data that are obvious in Figs. F-13 and F-15, this 
same frequency filtering process was applied to the 
first-stage edited data of all three segments of the 
Georgetown data. Figures G-9 through G-32 show the 
results. 


The "before" and “after” filtering times of 
arrival of the standard third-cycle zero crossing and 
the voltage at the 12th sample point were computed 
also. Their values, along with their standard 
deviations based on 8192 individual pulses, are 
presented in Tables 3 and 4, 


A comparison of the mean values of the 
second-stage edited data after filtering in Tables 3 
and 4 with the corresponding mean values computed for 
the unfiltered data indicates that both processes 
obtained the same result. We conclude, therefore, that 
the Georgetown problem is not the result of the 
coherent interference shown by the histograms. It now 
appears more likely that any measure of pulse arrival 
time is so sensitive a function of observed shape that 
even more sophisticated methods must be used to 
quantify Lt. Unfortunately, program plans did not 
permit further development of the analysis. 
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4. CONCLUDING REMARKS 


We hope that we have shed some light on _ the 


intricacies of Loran-C pulse propagation. However, 
much still remains to be understood. If there are 
others who wish to contribute their expertise, the 
recorded data can be made available to them. For 


copies of the data tapes, contact the Loran System 
Program Office of the U.S. Air Force Electronic Systems 
Division at Hanscom AFB, MA 01731. 
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Appendix A 


DATA PROCESSING COMPUTER PROGRAMS 


The fundamental data processing was accomplished 
in two stages, as shown in Fig. A-l. The first stage, 
referred to as the preprocessor, accepts the raw data 
tape as input, along with certain control parameters 
punched on IBM computer cards. The data on this tape 
were formatted by a Kennedy Model 9232 buffered 
formatter. The data from the fixed site were treated 
as the master system and were normally processed first, 
followed by field site data, which were time 
correlated, by program control, with the master system, 
Data fill of zeroes was output, when necessary, to flag 
missing simultaneous data at the beginning or end of 
the field site recording. Other than time correlation, 
both sets of data were processed in the same way. 


The output of the preprocessor includes (a) a 
9-track magnetic tape (condensed data) collated with 
all data from each mission, (b) a punched card deck of 
control parameters, and (c) a tabulation with means and 
standard deviations of the first 16 measurements of all 
50 pulse blocks and the total 20-min data run. 


A first-stage edit function is performed by 
removing all data that may have reached the limiting 
values of the A to D converter. These data are flagged 
and not included in any subsequent calculations. When 
it is found necessary to edit any measurement, the 
entire set of 32 readings of that pulse is eliminated. 


Finally, a plot of the absolute mean values is 
presented showing the approximate shape of the loran 
pulse envelope and behavior of the zero crossing 
points. 


The programs used for preprocessing the data are 
documented in Program Listings A-l through A-28 and 
Flow Charts A-1l through A-28. 


The second stage of processing starts with the 
second-stage editor. It accepts, as input, the output 
tape of the preprocessor and the punched card deck that 


contains the editing criteria. It then performs a 2 
sigma editing on each data point in the 50 pulse 
blocks. Again, the entire pulse is removed from the 


computation of new means and standard deviations if any 
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one measurement exceeds the threshold. Intermediate 
means of each page are printed at the bottom of each 
page, and at the end of each data set the final mean 


and standard deviation about the block mean are 
printed. 
The final page of printout contains certain 


group/phase measurement (GPM) computations. The first 
four values in the row labeled T10 through T1l6 describe 
the distance from the zero crossing to the measurement 
point for the 9th, llth, 13th, and 15th sample of the 
pulse. The last three values, labeled GPMI12 through 
GPM16, are intercept points of a straight line fit of 


successive five point sets of amplitude means 
surrounding the third cycle. ATL of the GPM 
calculations use the "page" means, A made up of 14 


, 
block averages. 2 


T(N) = - 3.18 x AY /(A, - AY) 


CPM(N) = 10 > la, [/c2clay] - la gl) + lal - a) 
i 
i = N-8,N-6,N-4,N-2,N 


The programs used for second-stage editing are 
documented in Program Listings A-29 through A-33 and 
Flow Charts A-29 through A-33. 


A secondary data processing system, the 
group/phase frequency filter (Fig. A-2), accepts, as 
input, the original Taw data tape and control 


parameters from punched cards. Although several options 
are available in the program, only histograms and 
frequency filtering of one pair of measurements were of 
interest. A specified amount of data is reduced, and 
the llth sample, which is near a zero crossing, and the 
12th sample of the pulse, which is near a peak, are 
presented both as a two-dimensional histogram (in 
pairs) and respective marginal histograms. Means and 
standard deviations are also computed. Both sets of 
data are then independently subjected to a Fast Fourier 
Transform. Any unwanted energy spectra are then notch 
filtered (coefficients set zero) and then returned to 
the time domain by use of the inverse transform. The 
reconstructed, filtered data are then redisplayed in 


= £6 = 


the initial histogram format, along with new means and 
standard deviations. 


The programs used for frequency filtering are 
documented in Program Listings A-34 through A-58 and 
} Flow Charts A-34 through A-58. 


| The three sets of programs described above call 
; : subroutines from a library. These subroutines are 
. documented in Program Listings A-59 through A-67 and 

Flow Charts A-59 through A-67. 
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Fig. A-2 Frequency Filter Data Flow 
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COMMON IMD, CROTUT,TOFER, TREEE, IQEC, TRC, TREE TTMR 
COYMON/GE TOS T/10°Ty TRUE 

COMMON/SITAS/ISLAS 


c 149020) LS US=O TO PICK 2NN C=T OF IA SAMPLES IN SUB BLOCK 
c Ist¢21) 1S USED FOR INTEOMED ATS OQINTOUT IN SUR BLOCK 

129=6 

12°F=0 

PEANCIROSLOGIDE ICH EK DK ade ® De CIHO CTD De Lal2elod, IBIAS e THO(20)— 

Vreoe2h) 

TECIHOCL2b D203 

2 THESBETREF (CL, LHMC13)) 


TSF eTREF( 
IsdCL2d=1T¥ 


THOCL3)) 
1, THDCL29) 


I+9(Le@)=IT 3eTHOC12)) 
THDCL3V=L TM 02, THICL2) 
fs 1 


62 10 1 - 
Fad 19S, L992) TREFB, TREE 


lang ECxMATOLLIS)D 
LON? FL2AMAT(2FL5.7) 
es. 


Program Listing A-3 Preprocessor Initialization 


INPUT LISTING AUTOFLOW CHART SST - LG2AN 1 CONUP PHASE 


(LIST, WaMsod 


ere CONTENTS 


SISZOUTINE HOCHK CICHKD 
OIMENSTICN ICHKCED » JCHKIL DS THOCS9) »PBUFI12e) 
COMMON IHD 
COYMON/GSETNAT/IPT 
DATA ICHK/* DAY", 
totes 
a? Tre cIPtT, 1603) 
Cc FAD HEACE® 82ND OFCODF 
CétL RHO 


rene 
YEAR ES IXMTA HS HSTTE TS EMODE NS 


oe JCHKITDD GE TM 1 
ASTTECTET S 1NOOVICHK( TDS THOITD 
eTITECIPT,ICOL) 

JOMK(LV=THOUTD 
1 CONTINUE 
RETURN 

10M0 E“IMAT(ZOK,AG,* OES NOT MATCH: 

IMOL FOIAMAT( 26K, MAILE USE CIDE ON TSP 

1993 foamati*ti) 

END 


*13,° TS ON Tape’) 
ih 


Program Listing A-4 Preprocessor Header Information Check 


INPUT LISTING AUTOFLOW CHART SET - LORAN 1 G2QUP PHASE POEPPOCFSSO® 


PreEPPOCEssne 


| 
| 


an 


«be 


V9/ 2B 474 INPUT LESTI’ AUTOPL OC @ CHaet 


FORTRAN MIT ULE CUIST.NAMSO) 
Caen Ny weer CONTENTS 
) cn 
2 FEL22 be tHOL SD) 
} SION SPOTOTO2, Led TORE TRE EF 
S OUTOT,T2EF9,7 ataee 
5 aTsper, Paue 
6 c 
, 
e c 
s c 
1 c 
rs c < 
1? c T AULA 
1 ¢€ Y #UL? 
te c 
BS € ae 
\- € BRaTiCN 
17 c Tree 
is i Tit 
16 ¢ tore 
2? c 
2 c 4 SOMME BETER MIDNIGHT 
de ( P TION 
> c 5 WTOUT ELAG 
26 c P 
37 C 
af c 22> CONT 
23 c =e SLACK COUNT 
26 
2) 3 
3? GO Ve (hy dele LLP bet 
ua ¢ TIMe TO EVEN SECOND 
36 & gePrecyy3) 
3 TFCTPER NEO) T9 4 
3* oo 1 
a7 LeKerl 
30 CALL SHFTL( ) 
an CALL PCOTAC IH TUL fy?) 
40 1 CONTINUE 
41 & CALL bATene 
“2 c aTVexvatl 
ca J=!4UF (11S) 
LS Call SHFTAE Lede 24) 
pa CALL SHFTI CL tee? 
ae Cah SCOTACIN SCS) ety?) 
47 c ChL(ozaTio T 
“e Je PoUFclis) 
50 COLL SMFTL OL Jee) 
oa TE(L.GT.O) 
<2 IhD(YL)=9 
33 ? Kel,7 
ae Ce SHFTLCL Jee} 
$s 2 MMVCEEESEmOCL Ebene 
S¢ JetSuF Clie) 
$7 CALL SHFTLET Jyh) 
$e CALL ECOTNGKs 192) 
«9 yelonewjyev 
4 c weae 
41 17) 
&2 e Ctbehe™D) 
&3 CALL SHETICIHO OAD Jee? 
he € wate cnur 
6s Call SHETLOLHM OAD deed 
ded € mone 
67 CALL SHETEULHD OSD ty 4) 
se 5 Jaber cui7y 
ge c SLT© cone 
7 CALL SHFTL OED e d+) 
71 CALL ECOTS(INS (a) 9 fe?) 
72 Cc weatHee COP 
13 CALL SHFTL (Ts de®) 
as CALL BCOTACIMN CISD Ee 2t 
75 € wore COO 
% CALL SHETLOIMOC IAD dee) 
7 € INSIDE Teme 
7e CALL SHFTL (Ted, 179 
79 c BCNTPLIPTM CIO) Te?) 
69 c Curstor reve en?¥) 
81 Jereueciies 
82 CALL SHFTLCT ote?) 
"3 CALL BCOTALIHOOM I hed 
4 CALL SHETLO Ty Jy 12) 
es CALl SCOTACIHICe de Ted 
Be Cc peessucr 
87 CALL SHFTLOIT Jet) 
#8 CALL PCOTE(IHDEED eT e?d 
£9 JetSuF (LL) 
pad CALL SHFTLG Te syt?) 
a1 CALL BCOTR (Ky E43) 
92 THC T=LHOCe TPL ON eK 
3 3 Get 
Nes CALL SHETEC Ted 5 1A) 
95 J=IBuF (120) 
9% CALL SHET(T Jy) 
“7 CALL ACOTSEIMOULTD Teed 
cr) ¢ 
9 . TEM® PATCH TO $sTaeT CN GMO VEFCI@D 
too ¢ TECTHOC2 De THD C17) NE6998) GO TO 3 
101 ¢ 
102 ¢ 
193 PHO CLE DETMO (D2 DOIB OME THOIT ID C6OeTHOUTS) 
19% @ETUON 
19s 106 weTTecrer om 
1ne ra to 3 
be? LL weTTECTPT,LOOLD 
19e $1oe 


“fe 


’ 


Gone owace 


pece 


7c essor 


Program Listing A-5 Preprocessor Header Information Decommutator 


RE eR KA ER 


52 = 


O97 28/76 


CARD NO 


109 
110 
tii 


09/28/76 


FORTRAN MCOULE 


CARD NO 


Omran 


09/28/76 


FORTRAN MODULE 


Carn ac 


ee 
mM ODRPYO PUN 


re 


-. 


ad Bite 


INPUT LISTING 


1000 FORMETEC® EDF ENCOUNTESEN, 


AUTCFLOM CHART SET - LORAN L GROUP PHASE PREPROCESSOR 


CONTENTS eens 


Tey AGAIN ///) 


1001 FORvaT(® NO OO CARO THIS FILE®) 


END 


INPUT LISTING 


(LEST Nase) 


SUSZCUTINE DATCOe 
Cc ay 
c ON fnToa FOUR SITS WERE 


OMM424 IMD 
N/GETOAT/IRT, SUF 
(114) 
115) 
SHETC LT ede) 
SH=FTCL Kyo? 
SHET( Lyd ye) 
TUF CLLS)ET 
T=et2J=(115) 
09 1 L=L5e120 
JETSUFCL) 
CALL SHET( Le dea) 
TSuF(L-Lb=! 
Cate (IyJ528) 
1 CONTINUE 
RETURN 
END 


Program Listing A-5 Preprocessor Header Information Decommutator (concluded) 


AUTOSLOW CHART SET - LORAN 1 GROUP PHASE PeEEPEOCFSSOS 


CONTENTS wees 


> Su3 TO COP7FECT FO? WIOING EVeens TN CONTROL PANSL 
INSERTED AFTER TI4® OF DAY 
OlmMeNSION IRUF(129),1THD( 40) 


Program Listing A-6 Preprocessor Data Format Converter 


INPUT LISTING 
(LIST NAYSCD 
ease 


SUBZCUTINE BCOTA(T ed K) 

M=0 

I=9 

Le-a96 

CALL SHFT(SeMyL) 

00 1 N=lek 

t=Te19 

CALL SHFTLC SoM e4) 
1 feley 

RETUSY 

END 


Program Listing A-7 


RR OE) ee wr One. oe) Ce 


AUTOFLOW CHAPT SET - LOPAN 1 GROUP PHASE PREPROCESSO® 


CONTENTS aeee 


Preprocessor BCD to Binary Converter 


26 


O97 28/76 


FORT? aN MODULE 


F caro Ne 


¥: 


RE RE OE RR TG 


(LIST,NAMSQ) 


eres CONTENTS 


SUBROUTINE PRNTTL 
Ipr=6 
wrrTeE (IPT, 2009) 
1 WRITEC LPT, 2002) 
WOTTECLOT, 160G) 
WILTEC IPT, L001) 
WRITECIPT, 1002) 
WOTTECIPT,109*) 
WITS CEPT, 1904) 
WEITECLPT 10059 
Werrecrer Ol) 
WATTECTPT,1996) 
WeITE( TPT, L007) 
WAETE CSOT, 000) 
WRITE (LTPT,1909) 
WOITECTPT £019) 
werrecrerT, lori) 
WOLTECLPT,2001) 
WaETE(TPT, 1012) 
WOLTE€IPT,1023) 
WOLTECIPT, 101s) 
WRITECTPT»{O15) 
WeITEiteT,iolAs 
W2ITECTPT,LOL7) 
WAITECIPT,2¢01) 
RETURN 
OOM FORMATO 32K TIS ae gE OH PHU Ea MIE OE MINN DIK S/H SKy 
L*AAAAAAT, EX, *PPEPPEO*,SK LLP) 
LOOL EORMAT (32g MSIE yey HHH HCH PHHY LAN, MULE oy MUU GBR GOSS 4g SKy 
LYM SK, AAT SK POOH LEH, COOH EY, MLL OD 
1OC2 FQOMAT (32% eS d 8 Gk g PHHMHHMHE 82K, FULT Sky SU THES 6K 
LSAAREABAE LK, PPPDDDIO® EK, MEL OY 
LOO3 FORMAT( 32K MSIE eae HME GK y PHT GaN MIU GOH MUI EX SLE TKe 
LY AA GK y FAAS eK GP PI>h TOM SILL OD 
OGG FORMATO DAK e SS SUP SIY She MHEG AK y THHE eh, MIE GM GMO She M/S Pe 
12x, "AAS SX AU GK PPO 8, TOK, LLY 
LOOS EQQMATEZTE pM SISISS 8g Shy PHM EG LPH ES Ky MUSIUUU g EKG F/G IX 
Leahy Gx, HAA OX, Poh oP OMe PLELLELEL 
1ODK EQ2MAT( 3K» "GGGGGG® EE. 12IZWOLP HAY» CCDOCTI" YS 
LEPPPPPPPH LOX eM //8 LU, HOPPE PDH SEX yHHHt ox, HY 
2ESSSSSS9 AER 
1007 FoaMaT(*® GE 
L8UU' SX _ PP 
2HRAt oe 
Looa FORMAT(* 
L'poPpePpppr sy 
2*sss5ss? 


PUL Sg FLIUY pa hy 
EK TARAAARY BKy 


€ Vy MV ax tI, cx, HOF 
EX,UPOU LOK OSS EEK I PEE aK MODE LK, PHY 
Ky PSEULEW, Shey, EEF) 

LOK, FPPDCR PAE, EW, HON LG, HOD SY 
gAV oO PPPPPPP! 5X y MHHEHHE HY yh Xy 


GK, UU oKy 
SK, HH ES oKy 


LEP eX FUL GK 
SARABAY SSX 


2tee* OK 
1010 FOeMATC 
LOM Ky MPPE GLI ESSE 2K yg MPO Ny LAX, MHHE GSK ye CHHY yoKy PAA GOK y AA Ny 
26K SSH yeXe S589 SX 
1011 FOeMAT(® GGcccc 


COI7DND § AM s MIUIUUU 9 SXo 


LOPP* 4 1LOX, *//% 99% ah MART OK, SAN LEX, 
2*SSSSSS*.S5Xo TECFEE 
L022 FORMAT (Gex, *TTITITITT 4X5 FEEV EK, PSSSSSSh SMe TTTITINT! D 
LOLS FQRMATIG7TY ITTY, T¥LPTe eit oTKetTTO) 


7%MW.sssss eA OTT 
econ, 7x, HTT AD 

LOX e "SS eSKe MSS oT OTT OD 
FECEE 05x, "SSSSSS*¢SXe TT") 


LOLS ECOAMATEST SITET, TY, 
LOLE FCOMATISIK, TTT? y 7K 
LOLS FEQSMAT(GTHy ITTY, IX, 
LOLT FORMAT(GT KS TT 9 ?¥e 
2000 FORMATCLHLy//) 
2001 FORMATI//) 
2002 FOOMATUS///) 

FND 


Program Listing A-8 Preprocessor Title Page Printer 


cae SA UT ha Te Rte BE tee = 


o 


SHPO LEST EN AUIDFLOW CHAST SET - LOREN 1 G2GUR PHESE PREPROCESSOR 


. 


4 


TE RR RRS 


09/28/76 
FORTRAN MONULF 


Caed NO 


09/26/74 
FORTRAN MCOULE 


Caen NO 


INPUT LISTING 


(LEST »NAMSOD 


1000 


INPUT LISTING 


CONTENTS 


SUBROUTINE DATE (JCHY) 
DIMENSION THUF(12°), THO(SO) 
DIMENSION Al ,B)eMOUL2) -ICHE EL) 
DOUALE PRECISION & 

COMMON [Hn 
COMMON/GETOAT/IRT, [9UF 


DATA MO/3L GSI eID LING VSL LID 212 e243 9273p ING 3345 ESS 


ret=s4 

TssCHRCLD 

TY=JCHK (2) 

00 1 J=lel2 

TFC E.GT.MO( JI) GO Tt 1 
IMJ 

THD (23)=1M 

1O=t 

IF(J.EQ.2P CO TO 4 
ID=1-MOt 5-2) 
THO(22)=10 

GO T0 2 

CONTINUE 

RETURN 

CALL OOTOA(4,IM,IN,IY) 
90 3 K=le6 
WITECTPT,1000N (ACK eID eSeLeeD 
RETURN 
FORMAT(LSX ye (4x,A°)) 
—ND 


AUTOFLOW CHART SET = LORAN 1 


GROUP PHASE PREPROCESSOR 


Program Listing A-9 Preprocessor Date Converter 


(LIST »NAMSO) 


sees 


CONTENTS 


SUB2OUTINE DOTOR(A,IM,ID,IY) 

OOUBLE PRECISION Ay@eC, 

DIMENSION ALS eOVe8(E611) C1540 ,0(12) 
EQUIVALENCE (COCLoLD,COLD De (CL e100 ,D01)) 


OATA C/* ONNMONO *,*An nor ,ton 00°, 
1000000 *,* 121 %%* 211 tet at 
2° LLLLLL',* 22222 #4822 22¢,¢ 22 
39222222224," 332233 4833 32",* 23 
4? 332333 rae £4 66th 64 

se Ga", ISSSESSSSH, ES +, "555555 
6* §55555 *,* G6AS446 Het "BABAK 
T® 666665 "y*77TTITTITEG? 77 
at 77 "," BERE2S *,vop esaae 
9° BBKaKB t/ 

DATA O/7* 999999 *,*#99 e9*.* g969999', 
1* 999999 *¢,° s/*,? 4/ ye tz 
2°47 " 

K=IM/10e1 

OC L T=be6 


ACT s1)=S(01.«) 
K=IM-LO@(K-l)eL 
Oc 2 I=het 

ACT e2)=S¢01,«) 
Oo 3 Telsé 


ACT ,6)=S¢1y11) 


09 & T=lyé 
ALT s4)=8CT ym) 
K=IO-10*(K-1)e1 
OC 5 f#1y6 

ACT SS)=B0T KD 
00 6 t=les 

ACT >TI=B81F eS) 
wKervel 

00 7 L=1yA 

ACT ,B)=801 ye) 
RETUEN 

END 


"00 00*,"00 
ASC; 1. 
os" 22 ‘ 

33% y* 33° 

SE OG GEEnAGE 

wg S58 

& ytb6 66° 
ek 7 . 

6 *,*e6 est 

U 994,499 
yt “wt ' 


AUTOFLOW CHART S*©T - LORAN 1 


GROUP PHASE PREPROCESSOR 


Program Listing A-10 Preprocessor Title Page Date Formatter 


= §§. = 


ame S eae O 


oe 


09/28 /7> INPUT LISTIN AUTOSLOW CHART SET - LOFAN Lt GROUP Pease PSEPEOCESSOR 
EMAIRas KOMULE (LIST NAHSOD 
Caen an seer ONTENTS 


n SuaenuTt PRT 
? DIMENSION ASETE( 3 ¢ 12) 9A OMT R269) AMONEC Se RD OTH (AOD, TSUFLZFD 
2 DIMENSION WINTCD) 
z DIMENSION WEATHIC 2) ©) WEATH? (2,10) 
s OIMENSION ANOTEL2,2) 
ie DOUSLE BRELESTON ACETE g@AUMT? AMINE pe RTHL WEATH? ASOTE 
7 CoMeoN HHO 
e CO¥ACN/GE T/ERT, LAH 
q DATA WIND/® LC aH | 
1 OATa wESTHIsS* CLEAR e,* CATTERE* SO TE SOS? 
i Leronen T#,t0 ong fy taveecesre T LDAT* SFOS SCUTAT 4, "TON BOOTS | 
re DATA WEATHD/#THUNTERSY, STOI4S SIN SHO") "HESS = 
. ret VENZLPILE e,efoecy RLLILE SS ESN OW pfs i 
x 1s 2 HALL **GRQUND © SStEOS COVE SER ‘ 
Ve aHSUNNY H e,tHaze | 
le Cava ANOTEs¢S LORS st 
1? Para aStTes# ‘ 
13 LPNs : 
19 : 
18 DANA THOt, tEE x 
a NANTUCKE' , °T Tae be 
oe © TAO“ sires o,f *% 
>. axxt JUPITER &yEMLET 
e ANTHC<E*9T ' ries wy 
26 CQAN ST®,MULATES %,°100 <4 C*,PESIUM > 
27 
> AMoDestparmagy & S¥yaave %,tSTUSY | 
22 *,'Cacracere ECUE', *PMENT TE 
TLLaz ty ty DATA | 
21 SKY MAYE : SPOR F BCHe ers 
2? sf 
33 [?T=4 | 
IaT=21 
WeITE (IPT ys 101i) | 
1¥=70+ITHOC2) 
| 7 W2ITECLPT e LODO) IH (230 e IHN Z20 6 1Y | 
| 2a T¥=1970¢1H9(2) | 
| 39 WEITECEOT,1GULIIHOC LD IY 
20 Jats eyed 
“y WETTE (TPT pLOC2N(ASITECT y J) of 
“2 J=THO(2 POL 
“ WILTECLOT » LOO3DCAKMTICT I) yg T=L 92) | 
“ec J=THO(S)e1 | 
se WOLTECTPT yp LOOSD(AMOMELT p Jb eT =1e3) 
“5 WOTTECIPT 1C12) THOS? 
“? HO=1HD (6) *.03386 
ee WRITE (IPT S1905)HO 
40 B2ITE( IPT» 1006) 
en rO="HO'T) 410 
re W2TTE (IPT, 160 7HO | 
ry HO=IHD(EI/10. 
52 W2ET= (IPT. 1003)HO 
©. W2ITECIPT»21009) 
55 TH#D(10)/10. 
Ss WIET= (12T¢1903)HD 
<7 WITTECTPTeLOLODIHDULID 
Se WITTE CT PT» LOLIDEHOCL2D «THOT 3) pTHOC LS) 
55 vel 
“0 K=IHOCLS)/19 
41 1F(<.GT24) M=2 
42 K=MIOIK,SDOL 
62 LEMONCIHOC LS) LOD]? 
6% Worr=e (ret, 101) 
< WIETEL LOT, LOLG) WERTHEIM KD) p 4152) 
AS WRITE CIPT,IOLTIIE THOON LD ee be?) 
6? WOTTECLPT,LOLS)WINOU™) 
42 WITTECIOTSLOLIIIHOULTD 
49 WALTE (IPT s1GLo) (ANOTE(Ns IHDC IS D6L) pNeLe2) | 
10 c 
1 c 
72 TECTHO(5).NE.O} GO TO 1 
3 wPETEC2L) THO 
L Je c 
75 c f 
rz 2 RETURN F 
77 1000 FO24AT(35X,"GOOUP/PHASE TEST PPEPSOCESSING 2K, 12y20°/% e120) 
7a 1DOL FO2MAT(FO4, 3OxKy PIULTEN DAY NOV fSy AK yh YEAR NOLO STS) 
[ 792 19002 FOPMAT( PO, 34Xky Ta PECOAHED ats, 2KN,? 
\ en LOND SOMATIC AOC, PTAANSMITTED LOCATED AT Paay 
a1 1005 FOPMAT( 10%, 24L,9NSTA CECNRAT MON 
22 1005 FOIMAT( 108,24, *ATHOSEHEOTS BPESSUOES,FTe le? MILLIBARS*) 
e3 190K FCOMAT( (08, 56%, FE XTESNAL TEWOT2STURE (DEGREES CONTIGOADED") 
a4 1007 FOUMAT(G5K,*W=T MULE, ES.1) 
pe 1NOB EQ2MAT(45¥, "ORY AILEY SES.) 
a6 LOAD FOAMAT(FO*, 3K, MINTEONAL TEMPEDATURE (CEGREES CENTIGRAOENNT 
27 1010 FOOMATNO%, 24%, 9P2ZED. TIME CALISE 4,119". OF -OF SECONDS") 
68 LOLL FORYATENL'S///7/) 
°9 YOL2 FOSMAT((O% 24K, FATTENIATION IN D8: 
» 90 1LOL3 FORMATE FO ae tSTAOT TIME ho IGe? 
ar 1OLG FOIMAT(@5xX,*CLOUN COVES: 
92 LOLS FOIMAT( 45K 486," GOCUND 
93 1916 FORYAT(//// 35K, te 08® ey 
* LOL? FCOMAT(G5K,"GENFIAL CONDITION: * 5248) 
P 95 TOUR FCRMAT( Oy tang THE ATHEOS TY 
: % LOLS FORMAT ('G*y3eKy (GENUS REPETETICN INTERVAL? "yf 4,*.0E-O6") 
97 —Nn 


Program Listing A-11_ Preprocessor Header Page Printer 


HA ie 


| - 56 - 


+ te 


SOT SS Dae heen ew B+ buen N % ‘ 


reas Soe Se eee Ss 


Ovs 28/76 INPUT LISTING AUTOFLOM CHART SET - LORAN 1 GPQUP PHASE PRFPONCESsOR 
FoSTRAN MODULE (LISTNAMSO) 
GanooNe ane CONTENTS eses 


L SUAQOUTINE MAINE 
2 OIMENSLON LPUF (12781 DAT (BON) THI(69) 
3 DIMENSION [T48(3,256) 
6 DIMENSION GROTOT(3,14) 
5 DOUBLE P®ECISLON TEFFI?,25) 
6 POUBLE PLECISION TA,T89, IN, GIt 
7 OOHSLE PRECISION GONTCT 
e DOUSLE PRECISION TOES, TReFE 
9 COMMON IHD GROTOT,TPIEES,TREEE, REE LAC, TREE VL IMR 
10 COMMON/GEINAT/IOT, raue 
i lets4 
12 1yF=0 
PY CALL CLEAZIGRNTOT,42) 
14 GRI=IHOC LTD ®NAL ON 
1s c IHD(37) 1S TOTAL NO. OF PF YECTED PULSES 
16 IHD (27)=0 
17 Cc IHD(38) ES TOTAL NC OF MISSING RECCROS 
re IHD (26)=9 
19 NVANTELS 
20 1pcT=0 
2. 1£2a=0 
22 Celt TIMEC TS) 
23 rec=1ecel 
24 TREE CLS IRCD=T8 
25 TMS Ly f8C HOUL2) 
26 LTR 62 e12C HO(13) 
; 27 1T42 (3, 1RCP=THOC LE) 
2e TSP=TE 
29 IFC T2EF.EQ.0) CO TS @ 
30 CALL FINDT(TREFS,TE>,GRTyTOCT, LIF) 
31 8 IFINPONT.LT.14)G0 TO 9 
32 WRITE( TPT, 1003) 
33 WRITE (TOT »IO06) 
2% NPANT=0 
35 9 CALL SETUDAT) 
36 10c T=0 
37 c 1HO0(39) IS USFO TO CONT MISSING @ECCROS 
ze 1HD(39)=0 
39 ¢ 1H0(40) ITS USED TO COUNY 84° PULSES IN A SLOCK 
40 THD (40)29 
41 To=TseelPCTsG?r 
42 TFCLSFeN=29) GO TO 146 
“3 10 CALL PLCCKIPAT,IDCT.P2CT LET) 
“6 TEC TOCT.LT.500) GO TO t1 
45 CALL STAT(TO,0AT) 
45 TREE (2,1RCV=TD 
47 THO(37TP=THOCSTISTHO( SO) 
“8 THO (2 8dD=THO( Fede THO( 79) 
49 TFE(THOC21).£9.9) GO TO 100 
50 TFA (GODTOT EA, LP OLE CE. SAN) FETHON 
£t 100 NPANT=NPRNToL 
32 IFUIREF.EQ.0) GO TO 101 
53 TFC OARS(TO-TREFFD.LT..OL) GC TS 22 
£6 IE(TO.GT.TREFE) GO TO 21 
$5 101 IFC TPCT.LT.7) GO Th B 
56 11 TS>=TSPe7"GRI 
ST? 1Pct=0 
sé IFC TER® NEO) GO TO 16 
$9 12 CALL GETREC 
4 60 GO TOCL3+L2417,18). fer 
61 c TEST FO? CHANGF NF PECOROING MCDE 
#2 13 CALL MODESSIMDSV) 
63 TEC ITHO(S)D.NE.MOSV) GO TO 19 
{ 64 CALL TIMECTAD 
65 Cc TEST FOR MISSING R°C025 
BA 146 IF((18-T87).GT..6) GD TO 16 
67 Cc TEST FOR INTERMEDTATE SFCNOn 
ae TF((TBP-T3).GT..901) GO TO 12 
69 TERR=0 
70 TaD=Te 
a} 15 TF(600-10CTI6 +8419 
7% 16 (FQe=4 
73 THO (39) =THOC29) «2 
4% Ty&=0 
7s GO To 15 
76 C 19 SHOULD ATTEMPT BACKSDECE ON TvOUT TF PCSSISLE 
7 19 CONTINUE 
78 17 TEC THO(S).NE.9) GO TO 2n 
ee TECT2EFLEQ.N) CO TO 2? 
; £0 T8=9999. 
e1 GO To 16 
F2 2L TE(TS.£9.9999,) GO In 2? 
83 CALL SKIPFL 
; 84 22 FNOFILE 21 
$5 20 RETURN 
a 18 WRITE(TPT, 1900) 
£7 stop 
& 1000 FO@MaAT(® NO OF Caen Ters FILES) 
fe L093 FORMATULHL LE (2% 6 MAVEL HSE) 
90 100% FQRMATUIK LAC IK SNEVED A/D 
ey END 
c 
z 
- 3 Program Listing A-12 Preprocessor Main Data Handler 


i 2 


% 


+ ate 


al Td ae Pa 
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4 arsier. Trew“ 
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to co Tor 
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12 KeLneCT-tbege fHO(2M) 
13 CALL HVCIOCT»<,04T VID 
is 1 NTINVE 
15 a 1H ( ts FMD PYTEOMENTAT= POTVTOVT Ce EACH BEADING 
Ye TE CLHO(2p).50.9) GA TO & 
17 me1OCT~1° 
16 acttetier WOPEDATIND eM Oe TOCTD 
Ye » IFLIDCY.* 50) GA 19 3 
a¢ 2 CONTINUE 
21 3 rectet 
22 eetue 
23 (NON FCIAMATCLSF S61 
24 ene 
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C279 N tere CONTENTS eeee 
1 ISPOUTINE ELNOTETEEEe TR, COT.TPCT SI JED 
2 IMPLICIT 2EsLeats-H“,r-2) 
3 MIMENSICON Tyr 0129) 
os COYMON/GE TIST/1RT, TEVe 
« THs ASS (TO>-T2-F8) 
7 4 TEUCTOWLT.~.92) PETUON 
7 OF 1 T=2e L0G 
I T=TO-Corel 
6 te (oaes(Td.LT..08) CO TO 2 
10 4 CONTINUE 
11 2 TecTevonth.T) 
12 NeI/7 
13 TFUTER-TREES) 3oS 510 
te 2 TECNLEQG.OL 2ETURIN 
15 oo © Teles 
16 6 CALL GETREC 
17 CALL Time treed 
19° © RETUIN 
19 10 TECTPCT.EO.9) GO TO tt 
re NeNeL 
t, a1 pect=7-yect 
22 ll Ta¥eTee-NeT#G2T 
23 1J= 
4% eeTuey 
25 €ND 
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sees CONTENTS eeee 


SUBROUTINE TIME(TR) 
DIMENSION LKUF(L2°) THC G4) 
MOUBLE PRECISION Te 
COMMON [HO 
COMMON/GET IAT /IRT, Late 
COMMON/ALAS/LEALAS 
MePSUF(Ln3t 
1820.600 
MOL Teled 
Jetell 
COLL SHFETIIL KO) 
CaLeL BCOTACMs Le?) 

1 THE (M-[HO( J) ) 86088 (3-1) 674 
CALL SHETLEL Kk eFd 
CALL FCOTS(Myt 2d 
TR=Teervs.01l0 
MSTSUELLIS) 
CALL SRETIOL KR ,20) 
CALL FCOTS( Mel 5) * 
TH=TH evs OOUOGOLNOCTATAS 
PETUON 
END 
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(LIST NAMES) 
eeee CONTENTS eese 


SHPEOHTENE SETEDET) 
OIMFNSION OFTOL) 
00 1 T=1,800 

1 DAT(TD= 9999. 
RETURN 
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SUBQOUTINE KY CTOCT Kk eDAToL? 
OIMENSION DAT(L),TAUE(1280, THM CO) 
COMMON IHD 
COMMON/GETOAT/I2T,129US 
T=1T8UF («) 

CALL THETIC Je 1y 3?) 

CALL SHET2(Jyly-14) 
OAT(INCT-1)=JeL.220°3756-3 
CALL SHETIOJ9 1,32) 

CALL SHET2( Jy e-16) 
CATCIOCTI=Je1 22963785 -2 
CALL EOITCIOCT,OAT SLD 
QETURN 

END 
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na 2 f=beeo% 
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GROTOT(3 eT DGENTOT (Ay TD eXYNGT) 
MYC TDEKYC DIZ CONE) 
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XYSDCTI=*SSITECVICTI “KV CT) e272) 
3 2 CONTIN 
23 WOTTE CLOT LOOT ITO SKYNOLD SEHD (SC Dy TR DE29) 
2s WALTECTPT,LIOIVIEY 
WOTTECLOT IOC LI EYSD 
TEC THOUS) o) 62 Th 3 
WEETEC2LPTO XY, XV SVeNYNS (MAT CL) L=15299) 
3 RETURN 
LOOM FOIMAT(IX, LAES 6) 
TOOL EQQMAT( 2K Lore seul) 
1003 FOOMAT(® ALccKe TIVE © "yFL2. 7% 
LLOXe "TOTAL rULSES KEPT Lary 
2IOXy NO. CF PULSES ENTTED at yl2y 
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END 
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SURRCUTINE CLE S200 WN) 
OPUSLE PRECISION MIND 
OO 1 LeleN 

1 X(T)s0.000 
RETURN 
FNOD 
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FOr T@ aN MODNILE CLEST,%a4S0) 
cacd NO oss CONTENTS tone 
1 SUAROUTINE EOIT(LOCT»MaT yd) 
c: OIMENSTON [© (122), THI(G0) 
3 OIVENSTON MATL) 
6 Covvon THO 
$s COMMON /GETOAT/IET, 18 
6 PECASS(DATIJOCT=-1)3.47.19.) GN Te 1 
7 TECARS(DAT(IOCT)) GT219-) GN TOL 
e RETURN t 
9 L TOCT=INCT-2045 
lo THI (6O)aTHI(ECI OL 
11 Oo 2 Jz1,8 
12 TCCT=1DCT+2 
13 ORTALINCT-1)=9 999, 
14 DaATINcti=s-99, 
15 2 CONTINUE 
14 OETYRN 
17 END 
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1? GODTOTE >, LV=DSLATE ONTNT(>,LI-GEDTCTEL Teer 
13 2 CONTIN 
ira T2G2OTOT(3,1)e65 
15 WOLTE CHATS LeGldd 
16 WAETE (LAT, LION EGIOTOT Leh ele leles 
7 WIULTECLAT SL OMLVEGIOINNEZ AD ebedyted 
16 WALTE (LPT, LOO DINO(37) 
19 WAITECT ET, 1CO-PIHOEIS) 
20 co tT 2 
21 
22 
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26 ADNsa> 7a 

25 WeIT= 17,1907) 40° 

26 WITTECT,LOOSDITEOTITE De Ede iA) = 
27 3 @ETUTN 

26 1M00 FOZ¥ETILKetofe ee) 

29 LOOL FORMZTEM*K SIF ened 

30 19N2 FORMATOC® SAN TOTALS CV=S%_Tbe* PULSES aee ss FOLLOWS") 

71 INO? FOAMATE//, 10%, "TOTAL NUMSTA OF CEJECTED PULS=S =*15) 

3? 1004 EOUMET (SS LOK, TOTAL qimveca MF MISSING OF ce7eT PECOROS =*,35) 
33 1005 FOPYE TILED e eg Sper Pee) 

ae €ND 
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SUBROUTINE MOTESSUMSV) 
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1 SURROQUTING PLOT 
? DIMENSION THOCS Oe PSEC LE) 
3 DIMENSION GROTOT( 3, tAD Mere ( 16) 
‘4 OOUSLE PRECISION GEOTFT 
5 COMMON THO ,GROTOT 
6 COMMON/GETNAT/ IPT, T8UF 
7 IrT=6 
e WeETEC TPT» 1000) 
J ayvrsoO 
10 90 1 1=2,l6 
1h TECAMKC LTA MASS(GRITOT( ER LIDE EMK=OLESIGONTOTLS TDD 
12 1 CONTINUF 
13 MX=AMXe1 
14 MN=55/MX 
15 OC 10 1=1955 
16 L=55-1 
17 IC=LO7TESS 7074 
le TE(MOO(LeMDD EOL) TCEIL YY) 86553641085 67046E4 
| 19 00 20 J=lsls 
20 V=DASS(GROTOT(LyJe2)) 
I 21 MARK ( J)=LOTSTH 1976 
| 22 
23 IF(I.NE“K) GO TO 20 
26 MARK (SPE P22T HA IEMOYOM 25M 251655 2-9 
25 20 CONTINUE 
| 26 1O WRITECTPAT,1NGL) TC »MArK 
27 WSITE (IPT, 1002) 
26 @ETURN 
79 1000 FORMATE ML LIX LEC OT e6Kd 
I 20 LOOL FQRYATC LIX AZ LK LEC AT AUD) 
31 1002 FORMATO L4Xs140"1' 96x) 
32 —ND 
Program Listing A-23 Preprocessor Plotter 
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09/28/76 INPUT LISTING AUTOFLOW CHAST SET - LOPAN 1 GOOD PHASE PPEPPOCESSOF 
FORTRAN MODULE (LIST »NAMSOD 
cared NO sere CONTENTS coon 
1 SURRAUTINE SKTPFL 
t 2 DIMENSION [BUF( 122) 
3 COMMON/GETOAT/TOT, TUF 
“ 1 CALL GETPEC 
5 GO TOCLeLe3e2)o18T 
6 2 WRITE (6,100) 
7 3 RETUEN 
e 1000 ED2MAT(LX) "NO ON CAOM THTE FILE (eK r>eL") 
9 END 
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FOR TOAN MOIULE (LEST NAMSO) 
Capo NO thee CONTENTS eee 

1 SUBROUTINE Ge¥ 
? TMELICET FAL eS (AM yO=2) 
a DIMENSION THD ESO), G29TOT(2, LAD. MUTOTD 
‘ COMMON fHO,eROTOT 
s Tetse 
é DO 1 Teles 
7 J=re(les 
€ L QUT CED =-CADTOTE Ly S-D DOV. ASE CRNTIT Le IP -GIOTITCL I-20) 
q OO 2 32253 
10 Jeles 
13 OuT(Js)d=0.0D0 
1? De 2 K=1,5 
13 Le2eClen) 
le 3 WTC IPEOUT CUD eOAAS(GONTITO DSL) 
5 K=2e(Iel) 
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20 TEC IHOC35).N=67) GO TT & 
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2? Go To $ 
23 & ALTE CPT. LonLyoT 
26 5 ex=¢ 
25 $x1=0 
ae T=-12.5 
27 97:19 Teles 
76 T=Te5 a 
29 SxS eOUT(DD 
30 10 SeT=SKT*OUT(T DOT \ 
31 AaSK/s 7 
32 M@=SKT/125 
2 Te-1265 
36 00 BL f=lys 
38 T=TeS 
26 LL CUTELD=9UT(T)-a- fT : 
27 WOETECEPTy 1002) (OUTET Ds E=294) i 
26 OETURN 
ag LDN EORMATESS Seeley? T SelPDe31GKe SOM POI2d) 
“0 wor 
bal 1092 HOM eH det DODD 
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caer NO eeee CONTENTS sess 
1 SURRCUTINE SHFTOL sy JeX) 
2 CALL SHETCAT  SeR old 
3 RETURN 
« END 
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FO? TIAN MMOULE (LIST NAMEN) 
Caro An sone CONTENTS Paes 
1 SUBROUTINE SHFTLIOL etek) 
2 1-0 
3 CALL SHETLU Ty deol) 
s 2FTURN 
5 enn 
| 
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Caen eees CONTENTS =ne | 
1 SUAIOUTING SHFETIUT Sek) 
2 CALL SHETA( Ty Sokol) 
3 PETUEN 
4 END { 
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CALL EOTT COT asf ONDee THUESH OTe TFLG) 
TECLELG.FO61) GO TO IE 

LENT SECA OL 

Tscetscel 

OO 10 yelel® 

TeweeNatacsels 

SuM(les UTC yest eTEMY 
SUML Cp eSUM EC We TEMO 

SUM (Ze TI ESUM OP a sD oTE HOE? 
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CONTINUE 
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LEDAT CLD Lat eA td 

TF (MOOMLONT east). FG,0) GO ToL 
GO 10 ? 

LF (NA,ED.01 Go TAY 

00 HOl felels 

SUMC DoT eSUM EVE TINA 

SUM (eet Sus Cae tI ZNS 
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CaRO NO CONTENTS 


109 210xe* TOTAL PULSES EDITED #*e13) 
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117 1009 FORMAT (FO eee t EDITING AVERAGES USED Be / 020K BF 9.50 / 625K ¢4F 9,5) 
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\7 1*NE te 'CaROLINA Lt Te} * REACH The 
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60 TF (KR.GT,4) “ep 
61 KaMOO (Ke Sr el 
62 LEMODITHOCITS) »1OdeL 
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6A c WETTE CLOTS LOLA) (ANOTE (Ne IMO(16) 01) eNOLe2) 
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3 c 
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7? 1002 FORMAT (HOt se 3euetOaTA RECORDED ATE s IK, 3AM) 
7A 1093 FoRmMarc(ro *TOANSMITTFR LOCATED ATEte Jee 2a8) 
” 100% FORMAT (FO *OATA RECORDING MODE 4s 3X6 3A8) 
60 1005 FOoRmary® TMOSPHERIC PRESSUREDHFT 16! MILL TAARS SD 
61 1006 FORmaT( "0 EXTEMNAL TEMPERATURE (DEGREES CENTIG@A0E) *) 
62 1007 FORMAT (45x—*wrT HULA OFS,I) 
63 1008 FORMAT (4See*DOY AULA SFS,1) 
as 1009 FORMATCHOte eueFINTEMNAL TEMPERATURE (DEGREES CENTI GRADE) ty 
85 1OL0 FORMAT (HOt e Ione PREN. TIME CALTAI e LI Le? 06-09 SECONDS") 
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67 LORL FORMAT (HL e7s7y ‘ 
8A LOL2 FORmMaTi +O deus ATTENUATION IN DAEYEDS) 
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90 1OLS FoRMaT(4s CLOUD COVERT ts 
91 JO1S FORMAT (454 648," GROUND ’ 
92 CLOLG FORMATIZ/7/70395 Rte Oe 2Ahy  Peeeeeeeseet) 
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SUBROUTINE EDTT(OATarkYeTHeleIFLG) 
IMPLICIT REAL CS (A-C#b-HeO~Z) 
OIMENSION X¥ (16) eDATA(16650) 
IFLGs0 
00 1 yelel6 
TF (DABS (DATA ( je T)@X¥(J)) eGT.TH) GO TO 2 
1 CONTIAUE 
RETURN 
2 IFLG#l 
10 00 3 yalelé 
1 3 OATA(J01) #9999, 
12 RETURN 
13 END 
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RETURN 
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SUBROUTINE Gur SM OUT eRNT) 
IMPLICIT SEAL OR (bane OmZ) 
OIMENSTON SMe TAD eOUTIT #100) 
00 1 felee 
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‘ rer / 
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e@eeeteeeeeeee 
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1 Bec ewee eee seeee ee eee | 
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eee ete eee wee nwee oeeeeeeeeaee 
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eoeeeeeeeeee 
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oeee-® END OF DO e 
e Loop? e 
e . 
e e 


03,20--->} 
1 ' 2 
Were ween seceeee 
‘ 4 
4 wRITE TO DFY # 
, rT 4 
4 WIA FORMAT ‘ 
/ 1001 4 


Flow Chart A-29 Editor Driver 


= 102 = 


o9s27/76 AUTOFLOW CHART SET = LORAN 2 GROUP PHASE EDITOR PaGE 62 


CHART TITLE = PROCEOURES 


4 2412 Oe nnn ne wen nnn een n nent 
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secre rer tett eres ' * 12 t * 20 ' 
/ 4 [ e e ' e e 4 
7 READ FROM DEV / ' . . t ° . 
‘ TRO ‘ ' TRUE © . 1 . © TRUE 
7 INTO THE LIST / ' -* IFLG .fo. 1 © ' © TCNT ,€0, Q @rn--n0 
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eae e were cence es ‘ ! ° . ' e . 1 
' ' ' e ¢ ' ie 8 t 4 
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1 NOTE 02 ' ' IFALSE ' 1FALSE ' 
eeeeeeoeeveee 1 1 ' 1 1 ' 
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eereevreer eee e 1 1 1 13 1 1 1 ; 
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e e ' 1 Oa eww enn nen e rere n ent 1 . eee weer ene 1 
° ° ' ' ' ' ' t q 
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. ° f] 900 ' 10 ' eeeeveveeve eee ' eeeeveeveveeeeve 1 7 
. . ' weet te eeete reer 1 © 6FGIN PO LONP ® 1 © REGIN OO LOOP €8 I 
€. & ' 4 walte 10 OFV 7 ( © 10 58 te 16 . } © 200 © dp U6 5 - 
. 1 7] yer 4 fl eeoeoeveeeee ee 1 eeevreeeeeee { 7 
1 NO ! 4 vie Formats ' ' ' ' ' : 
' ' ‘ 1006 7] ' [depecceecscece 1 O2.25--=>1 ' fr 
' 1 7 FaoM THE LEST # ' ' 1S ot ' ' t 
1 ' eter went ereree 1 Ow eww wren rene nen n en n® 1 ' e-- onee f HH 
' ' ' 1 § TEMP © OATACyet) ot ' ' ' SuM(Lel) © ' ' ' 
' ' ' ro ' ' ' 1) SUM(LeEIZTONT ( ' f 
' ' | NOTE 21 (a | SuMilegd ® t ' ' ' ' t 4 
° 06 ' OFS 8S S18) 6 i) 1 SuM(leg) « TEMP of 1 1 i) SUM(2eT) ® i} ' * 
eres oe List «© To ° r * t 1 ot Sum(aetastenr ot ot ; 
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09/28/ Te INPUT LISTING AUTOFLOW CHART SET =~ (IRAN 3) GAOUP PHASE FILTER 


FORTRAN MODULE (LIST, NAMSO?D 
4 CaRO NO eeee CONTENTS eeee 
a 
1 IMPLICIT REAL®B(A-C sE-H 0-2) ; 
2 OIMENSIN 8125642) ,AVGP MP (25 ; 
3 DIMENSION JCHK(S) 10160) of BUF CL2AD - 
* DIMENSION GROTITC 3,320, AR P( 619452) 
s OIMENSICN OCE(8192) : 
6 COMIN IHO,GROTOT, ARR OQI AB 
7 COMMONS GETCAT/ L&T, 1 BUF 
8 COMMONS CARD/TTAsSDTOASNYFS 4 
q COMMON/PRINT/IPT 4 
10 1PT=6 | 
ut trRcs0 
12 TFLG#0 
13 TP =62.5 ¥ 
16 PHI=0.0 
15 €CO=2.0 A 
16 1b CALL CROITACICHK,AEFD § 
uv - TEC SCAR (19.69.999) GI TC & 
16 Catt INT 
19 Ttry=0 
20 2 CALL 290 
2. Cart CNCFCITRYD 
22 TFUITAY.EC.9) GO TOR 
23 TECHET?vL&.T) GO TO 2 L 
2 ENCEILE 21 
25 Cac SH ieee 
26 co rot 
27 6 CALL POATTCEIHI,LOT) 
2a CALL PR NT HD 
29 CALL MAINE 
30 Cace CONF ETAOU) 
uu Cau srati > 
32 S$ Care “Ist 
33 Care PLOT 
6 PFC ICR I20.69.09 GN TO Y 
35 TFC sCreI9d.FU.Od GO TO 6 
% cauc Toa 
” CALL ave 
36 $ CALL OF FT AC ARH, IME 25), 0) 
9 CALL OF FT ACAMAC LS EET E25D Od 
40 CALL PLOT ICA VSP We) 
a TEC ICH E9D.C.0) GO TO 7 
42 CALL EDITFC AVG ae) 
a) CALL PLITICAVIO WR) 
“ TO CALE PFFT ACARI, THOE259, 2D 
5 CALL OFETACRIRE LAD THE OD 92d 
46 THC SC EID.FQ.0I GA TO 3 
7 at Caue MIST 
an x Catt “vane 
ag Dd CALL PENTTCOTHD OTD 
= so CALL PRAT 
ais st CALL SNES CICHE AFF OED 
” $2 $ Cauy Toa 
$y Care axe 
s4 TECIC CMD HE ON OTTO LE 
ss CALL SUFS CICK AFF AOD 
56 METTLER LOD IDIMIUL DS TITAS STON: SES 
s? WEITECT A CO2D TP PHT ECO .T ADS 
38 co rot 
39 © MF w IND 22 
60 WOTTECTstbcctd 
6. stop 
62 1000 FORMAT 16,671K.9) H 
e3 LOOL FIRMATIIK 6990") j 
os 1002 FO2MATIGFIS.5) : 
es €ND ; 
' 
: 
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09/28/16 Tee uT LISTING AUTOFLOW CHART SET - LORAN 3 G®OUP PHASE FILTER 


FORTRAN MOOLLE (LEST, NAMSOD 
caro NO CONTENTS ecoe 
1 SUSPCUTINE C@OLNCICHK, AEF) 
q 2 TMPLILET REAL#4(A-C ,E-4 0-2) 
4 3 DIMENSTON JCH<KO05), [8UF(0 128), 1HOO SO) 
* COMRON IHD 
5 COM#IN/GE TOAT/I AT. 1 BUF 
6 COMPON/ PRINT/ IPT 
7 c L = POWER C©& Tw? (999 STOP) 
s c 2 = FIUPTER CL VES,O#N)DD 
9 c 3 = FILTER LEvVELENOT FQ ZEeC-aAwP CALC AT T12) 
10 = @ = T 2EACl FOR AMOLI TOE CALCULATION 
uw c THOC29) «= SS OAY ““umBER . 
12 c AEF © ANT EFFSCTIVE HEIGHT 
13 1PO*5 
itd B REACLIAD, LODINCICHK (LEDs Lede SPy LHIC 29) oAEF 
is TEC ICAK CLD LE 2.9) SCHRELD 12 
16 TECICHR (ODES .O) SCHKIGD AOL 
uw TEC ICH US DESI SCHKAIS HLS 
1s TFCAE=.EQ2.03 BEFe. 007 
19 THOC 25) =sCm 15) 
20 THOC2T) waCHe etd 
2t PHO( 29) e200 mm CLD 
22 TF CICH OS DLNE LOD THOC26 9812 
23 L RETURN 
24 1000 ATLELEFO.2) 
2 
Program Listing A-35 Frequency Filter Card Input 
9/28/76 tNeul LISTING AUTOFLOW CHART SET - LOR4N 3 GROUP PHASE FILTER 
FIATRAN MODULE (LEST,NAMSUD 
caro NO eee CONTENTS eeoe 
Ld 
1 SUReOUTIN aur 
2 PMOLELTET PEALeACa-H, 2d 
3 DIMENSECN EMCES DD, GROTITE 3, 220 
4 COMMON THO.GRCT OT 
s THO251*0 
6 c teee3s ULSE COUNT 
? THOS) « 
8 Cart CLeaetcerrsr,i92t 
9 c THOEVYD ES TCTAL NO, OF REJFCTEO PULSES 
to THOC STD #O 
u ¢ THOU33) FS TOTAL NO OF MESSING BEC EDS 
12 1HOC34D #0 
LB) RETURN 
“ ENO 
Program Listing A-36 Frequency Filter Initializer 
i - 112 - 
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09/28/16 INPUT LISTING AUTOFLOW CHART SET - LIBAN 3 GROUP PHASE FILTER 


FORTRAN MODULE tL IST,NAMSOD | 
Caro NO eece CONTENTS seer 
} 1 SUBROUTINE ROHO 
2 DIMENSION IBUFL L128) e1HO1G0) 
3 OOUBLE PRECISION GROTOT (3,32) 
of COmMMIN THO,CROTOT | 
s COmMMOW/ CETOAT/IRT, [BUF | 
6 COMMON/ PRINTS LPT ] 
7 : L=JULTAN O4Y 
8 c 2=VEAaR-19 70 
9 (- 3=XKMTR NO. 
10 $ @=SITE NI. 
13% c S=Cata “OCE 
12 - O=ATMOS PRESS 
3 c T=OUTSICE wET BULB 
le c B=CUTSIDE ORY BULB 
is ¢ PAATTENUA TION 
16 c LO=INSIOE O8Y BIAB 
wv c LLsOBJECT OF RUN 
18 : L2=HIU2S5 START TIME 
19 c L3=MINUTES START TIME 
20 4 IS=SECC WS START TIME 
2k c LS=wEarreR CODE 
22 C LO=NOTE CCE 
23 € LT=cat 
24 c LB=STazT TIME IN SECIS AFTER MIONIGHT 
25 c 19=4O0E OF RUN 
26 c 20=41SSION OF RUN 
27 c ZL=RUN NIMHER 
28 c 2250AY NUMHER 
29 € Z23eMONTH NUMBER 
30 c 25eFFT POINT FLAGISUA AY) 
at € 26=44P CALC CN 3RO CYCLE UP PEAK OE PULSE--=SFT tN CRDIN 
32 c ZESFILTER LEVEL MULTIOL TER USED IN FOL TR-----SET IN COPLN 
33 (4 2eePOWER CE 2 FOR NC. CF PCINTS TO PXOCESS 
vw c Z9ePIDIESSING DAV NUMA R 
Bhd c JOePULS® NIMAFO 
4 c JLsPULSE COME FABLE STaer 
” c 35ePU_SE COUNT 
ua c 3MMISSING PMECORPN COUNT 
39 c 40=58) PULSE PER ALOCK COUNT 
40 c Se IP FIAST BECIRO 
at 2 CALL GE TREC 
42 GO TOt se2ol2e900) -18T 
4) 3 CALL GCETREC 
0 GO TNH oD ehOoth de IRE 
4s c TIME 12 EVEN SECOND 
a6 © CALL ACOFAULL Se Le Pe MMUF L2G e THIOL LD 
“7 CALL ALOTAELL I, de2 eTAUE FIM oEMOELIDD 
44 CALL BCOTPELL Ss Se2e AUF pL 2MeIMNC LSD) 
49 c ATTENUATION 
7 so © CALL BTOTPELLS Ae 2 eT AUF LZR GIMNC DDD 
7 SL c OnJFCr 
sz CALL KCOTEEELSs 34625 TAUF LZR oEHOULEDD 
5) c MONE 
54 CALL ACOPPELLSs 5Se2e TRUE eo L2H eTMOUL IVD 
55 ¢ MESSEON 
se CALE RCOTBULLSs Pe 2elAiF ol 2% otMOL20N) 
s? ¢ RUN NUMBER 
sa CALL AIOTSCL Le by boi Mtie Gh2% eIMOCZLDD 
59 c PULSE NUMRE® 
60 CALL HCOTHELLGe2e2 efAUE L2AeTHNESIOD 
ot c JULIAN O48 Y AND YEA 
=~ 62 CALL ACOTAELL 654, 3e MUFF SL 2A eTHOELDD 
63 TECTHOU2Z9D.EQ.9) EHOC29 de TH CEL) 
64 CALL ACCTALLAGe Tye hel dtm LAA etHOCII 
es c xure COCE 
&6 CALL ACOTPU LLG s Sy Le EILF oLZBeEMOUSOD 
67 c “Q0E 
66 CALL BCOTPELLT ele le LW 2A eo IMOUSDD 
69 c StTe COE 
To CALL ACOTECLAT 2 2o CAUE LAV EMO SES 
m1 ¢ WEATHER CCOE 
12 CALL SCOTSHCLL I 4, 2efAUF oL 2H eIMOCLSDD 
13 c MOTE COOE 
” CALL BOOTRELE Te Syl ofIUS LIA eIMDCLODD 
% Cc INS TOE TEP 
1% CMLL ACOTAELAT S Ped e TAF L2H eIMOULIDD 
’? (4 OUTSIDE TEMP CWE TOY) 
Li) CALL BCOTAILL 2, Be IMU LAA GEHOET OD 
9 CALL BCOTECLER SS, 3eLAUF ph 2A oIHDCSED 
ao Cc PIESSUIE 
- aL CALL BCOTPULLA 3. SefAUF L284 THOT ODD 
82 . Gar 
63 CALL ACOTAULL9s 5, Se ISUF E26 eEMOUL TOD 
a4 c 
es c 
as c 
a? THOULADSERCELZ DS IHOOFTHOL LI DeSOeEMIEL SD 
es RETURN 
hy, ag 12 werrecrer,.to00) 
$0 60 TO 3 
1 At warrecrer, toot 
ou 92 stoe 
$3 12 wRITeCcrPrT.too2) 
" 62 719 2 
5 900 weTTEC TPT,» 2000) 
96 sToe 
cia 1000 FORMATI® EDF ENCOUNTERED, TRY AGAING o///) 
§ 8 1001 FORMAT(*® ND OO CAPD THIS FILET) 
99 1002 FJ TC? EOF FLAST REC Sutert 
\ 100 2003 FORMAT(C*® NO OO Ca®20 FOR FI@ST FILE SRIPTD 
4 lot ENO 
i Program Listing A-37 Frequency Filter Header Information Decommutator 
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nh AOR NR GABOR Tred OPE 


la 


09/28/76 


FORTRAN MOOULE 


Caro NO 


Cer~eOveune 


09/2 4/ To 


FOITRAN MOTULF 


caro NO 


Cant vruen= 


INPUT LISTING 


AUTOPLOw CHART SET - LORMAN 3 GROUP PHASE FILTER 


(LEST»NAMSQD 


10 


22 


iB) 
10090 
1092 


CONTENTS eeee 


SUBROUTINE COOECI TRY) 

DIMENS TOW [HOC GO), LTAB( Te SP eNTABITD 

OATA [TAB/Lelelelelelelelelerbe-bebebe-bebe~belerbebertele 
Um Derbe~ berber be-be-ts 

COMPCN IHC 

COMMONS PRINT/IPT 

Call GETTAAINTAB) 

00 4 I#1,4 

OF © Jehee 

O00 2 K#J5,7 

Lsa-Jsel 

TFENTABILI.NESLTAR(K, TDP GO TO 4% 

CONTINUE 

00 3 K*1ly4 

Masee-b 

THO(Re3O) IT ABIM, LD 

60 To 10 

CONTINUE 

WRITECEOT 1009) 

TTe vel tTRvet 

RETURN 

TF ¢THOCS) .6T.2) GO TO 20 
TFOMDICTADE 301, 224660.0. ANOe T6EQ.31 GO TO LY 
PECL cE Qed ANC ST HOUIOD SF Deke CO THE SOD.FO6S) GO TO LY 
TE CTS EQ ee ANOWEHOUIO). «CHC THOCSIDCEQLTI GI TO 13 
WRITECTPT LOO2ZVCLHOUND, Lede) 

1vev=) 

RETURN 

TE EMIICTHOC 300, 2066920. ANDLTLEQ62) 50 TO 1S 

TECT FUL DePHOCSOD ECL eH IMOCIONLFQ.9P GI TO LD 
TECT EC. DePHOC IONE Qe STH CIMT NOPE SLED E GN TO LD 
PEELE Qed ANDO MINI TROL SOD 2DFQLLI GO TO LS 
WOTTECTOT SLOOZD CEHOCND  NaSL odod 

tte yeo 

RETURN 

FOAMATONL 20K, "PHASE CODE ROT TN TABLES) 
FYOMATE PL 20K, ENC ORPECT PHASE CIDF ENPUT* 6416) 
END 


Program Listing A-38 Frequency Filter Pulse Phase Code Decoder 


Ivrur CESTING AUTOFLOW CHART SET = LORAN 3 GROUP PHASE FILTER 


(LEST, NAMSOD 


CONTENTS @ree 


SUPSCITINE GETT AR(NTARD 
OPMENSTON NT ABC TD SIT RDCGOD SE PUFUEZAD 
cco THO 
COPMON/GETOAT/ 1 OT, TRUE 
COMMONS PRINT TPT 

09 Lb tle? 
Jeloett-Lbelo 
Letsuety 

CALL SHFTLIK st 4854) 
CALL SHFT ACK ele -LOe4t 
NTASCTD eR TAS SORE 
RETURN 

ENO 


Program Listing A-39 Frequency Filter Phase Code Generator 
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8 ENR aig TAS RT TR TEER Be toe _ 


ee 


a AT IT 


09/28/ To INPUT LISTING AUTOFLOW CHART SET - (IRAN 3 G&2OUF PHASE FILTER 


I 
\ 
FORTRAN MODULE (LIST, NAMSQD 
Card NO ones CONTENTS 
‘ SUBROUTINE PRNTHD 
f 2 OLMENS TON 4006 (3,8 )¢ THOU 40), [9UFC 128) 
I 3 OLMENSION WINDS 2) 
4 OIMENSEON WEATHLO2—5) WERTH 202,10) 
; s OIMENSION ANOTE(2,2) 
| 6 ODOUBL= PZECISION AMCOE eMEATHL pwEATH2» ANOTE TS 
T DATA WINO/* LOW', THIGH TS 
8 OATA wEATHI/* CLEAR *,"SCATTEZ=*, "9 T) 5SOR*, 
q L°AROKEN T*,90 90% Se OVERCAST®4® TO LOOk* ,*CBSCURAT® »* LON 199847 


10 DATA AF ATHI/*THUNDERS TORS e*RAIN SHI 
aL . UORTLZILE *o* * s*FREEZING',® ORIZZLE* 


wERS ts 
SNCw bd, tod bi 


12 2*hate "GROUND FY,*0G ©,*FOG COVE® te "e 
13 3* SUNNY ’ M4 € 7 
4 DATA ANOTEsS*SES UNG® "NC PRCBL®.PENSee 
is OATA AMCOEs* %* *s 
16 L®LORAN C *e*mASTER . 
tT 2°LO84aN C SECWWCAR® te 
18 3*LO2AAN O MASTER ¢ % 
I 19 @*LORAN O x OR 2 St e*LAVE es 
20 S*LNPAN O t Stave fe" A 4 
| 21 COMMON IHC 
22 COMMINSGE TOAT/IRT LAE 
23 CIMOWs PRINTS IOT 
24 fwT=21 
25 WETTEC LOT tOLle 
| 26 = Iw=7TQeT OC) 
| 27 WELTEC LOT LOODDTHNC IDs THOU 229 EY 
2A TALL DATECLHOC 29D TVET) 
29 WEITECLOT 19200, IN, IY 
30 Ty=l9roetHoted 
gt WPETEC TOT Loot Tote Y 
32 MPETECTOT L902) EHO(ed 
3 weLTECtPT,LIOs Loy 
% de THOUS Dod 
ch) WRTTECLT OT LOOCD (AM ECL eS Do Tele dD 
36 WETTEC TOT es LO22PIHIE SDD 
” WETTEC TOR sLOL2ITHNGGD 
aL) HOSTHDCAD&.O33A4 
39 WOITEC LPT, IOSD 
40 wWRITEC TOT, L004) 
at MO eTHOC TD /LO. 
42 WEITEC TOT, LI0TIwO 
4) HO =IHDCAD S106 
as WEETELEOT «SL IOADHO 
“5 WETTECL OS 61009) 
46 HO aTHICLODZLO. 
4? WOLTECLOT, LDR DED 
4A WETTEC LON SLOLOD EMC ELL Ds PHOCLIDS PHO O20, THO 2Ld 
a9 WRLTE CLOT SLOLUL LAV ULE CGR RE SU MCRSS 
$9 Mel 
st KeTHOCLSI710 
$2 TECKLGT OT Me? 
$3 Re OOCH SOL 
4 LeMODEEMOCESDLODOL 
$5 WOTTEC LOT, LOLAD 
$4 WOLTECL OT S LOLSD (we ATHLON RD ede 2d 
st WRLTELT PT SLOPE UWE AT CEN LD eb edd 
"8 WRTTEC TOPs OLS at Ee) 
™~ 39 wEITECL PT LOLVd Mott) 
60 WETETECLPT ALOLAD CAN TEENS THO CLODOL De Nebo 2D 
61 TECEWOO3SD.CT.ON GD Te 1 
e2 id 
€3 ¢ 
64 WPITECZ UPTO 
45 CAcl TIME (TH) 
66 WRTTECEPT » 1021) TH 
6? c 
66 c 
69 1 Re TUaN 
70 L000 FORMAT (G5, *LOBAN SEGHAL ENWISONMENT® p Qe eh2e2Q PF el 200 
7 LOOL FORMAT PO oar etOAY Net ol AK et VERR NOOO ISD 
72 L002 FIeMaT('D 4, "OMTA CILLEC TIM SETE CCOE at gett 
13 1003 FORMATO ND 66x, *TEANSHETTER CONELGUOATION CIDE #%e 12) 
™ 100% FORMATO PO SSK, MOATA POCOEDING MDE  o3ee IAS) 
a hed 1005 FIAMATE 9D eo G4, PAT SPE RIS OSESSIPES "FTL e® MELLES ACS 
1% 1006 FURMAT( SD * 4Or,"EXTEANAL TEMPERATIRE (CESIEES CENTIGSACE DTD 
7 LOOT FORMATE S5 Ae WET Bi At e® Sold 
" 18 100% FOPMAT( 55x, *O2¥ PULSE SeLD 
9 1009 FORMAT(NOT eGR, * INTERNAL TEROSRATIIFE ECEGESES CENTIGRANEDYD 
80 LOAD FIRMATE SOT SOK, SONIELT TIM DE AML SSTINZRUN® ooo SUL 2 e872 De TLD 
a LOLL FOAMATENL' //) 
a2 LOL FOIMATINO® Sen, *ATTENUSTION EN O42 F E90 
a3 LOLS FOQMATEC TD", 44K, *STAOT TIMES 6G? ehrtet2y : 
a4 LOLS FORMATEISSX.*CLOID COVE? 
65 LOLS FORMAT(55%,46," G2 DLNU 
“> 86 LOLS FORMAT(I// 45k, eeneeeeeoe tesee sn 
a7 LOL? FORMAT (SS K,"GENEVAL CONDITICN: *,288) 
, BA LOLA FORMATE FOSS 4K, PWEATHED IED 
; 89 LOLA FORMAT('O* 64x, "GOP PEPETITICN [WFERVAL: Po 14, %.0F-068) 
4 90 WO20 FORMATE PO SOK PROCESSING CATE 2H b2e2C OF tg fZdD 
1 LO2L FORMATO 44K, TIME OF FEAST ALICK IN SEC S2* FLOM) 
92 VO22 FORMATO NO 4OK, "PULSE NUMMERD OS TS) j 
93 END 


Program Listing A-40 Frequency Filter Header Page Printer iq 
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09/28/ To INPUT LISTING AUTOFLOW CHART SET - LORAN 3 G@QUP PHASE FILTER 
FORTRAN MODULE (LEST,NAMSQD 
CARD NO ees CONTENTS 
1 SUBROUTINE MAINE 
2 OIMENSTCN [BUF L286) pCATILO001,1HIIGO) 
3 OIMENSION GROTOT( 3,32? 
4 DOURLE PRECISION TB, TBA +TO. GI 
s OOUBLE PRECISION GROTOT 
6 COMMON IHO,GROTOT 
? COMMCN/ GE TOAT/ LRT, LeUF 
8 COMMON/PRINT/IPT 
9 GR L=tHO1L7)*.000100 
10 tecrso 
i 1ERR=O 
12 Catt Timecred 
13 TeP=TB 
14 9 CALL SET(OATD 
15 tocT=0 
16 ¢ 1HO(39) 1S USED TO COUNT MISSING RECORDS 
7 140139) #0 
16 c THO(60) 1S USED TO COUNT BAD PULSES IN A BLOCK 
19 140140) 20 
20 TO=TaP+lPCTeGat 
21 10 CALL ALOCKIDAT, LOCT, TPCT TERRD 
22 TECLOCT.LT.1609) GO TO It 
23 CALL WTPUTITC.CAT) 
24 [HOCSTD SEHCCI TE OTHOIGO) 
25 THO (38) @LHOC3AD OLHOL 39) 
26 TF CTHO025 0.68 .EHO(280) GD TO 30 
2? LOL TFCLPCT AL T.7) GO TO 9 
24 ut 
29 
30 LFCTE2®@ NELI) GO TO 16 
ry 12 CALL GErREC 
32 GO TOLL Ve l2el Tol Bd, LAT 
3 ¢ TEST =)R GF OF RECORDING Mi IE 
6 19 CALL ACOTACLETs Le be CAUF eh 2BeMNSVE 
35 TF (THO(S) .RE.MOSV) GD TO LT 
» CALL TIME (TAD 
. c TEST FOR MISSING RECIPD 
38 le LFCCTA-TAPD.GT.AeGRTD GC TO 16 
39 [-! TEST FIR INTEQMEDLATE RECOWO 
40 TFCCTAP-18).GT..001) GO TO 12 
4. TERR ED 
42 TAPeTB 
a 15 LF 16)0-10CT1949,10 
“ 16 [Feral 
“as THOCS9DHTEOCSID OL 
46 WRT TECTPT tool) TAP 
4” on 13 15 
~ 4a AT wRETeC ter ,too29T0 
“0y fel) 
$0 23 TECTHDE 257.68 2001) GO TO 24 
sh tet-t 
92 TFECT.E 260) RETURN 
” GO To 23 
ry 26 [HOC25) =f 
35 ENCFILE 21 
So 20 RETURN 
7 30 CALL SKIPFL 
ra 32 3 tPF 
“9 18 wALTEC TPT .1000) 
60 stop 
6. 1000 FORMAT(® NO OM CA20 THES FILES) 
62 LOOL FORMATE* RECORD MISSING AT TIME* SFIS.) 
@ O02 FORMATE FO* 44x, TIME CF LAST AL ICH EN SFCSE* oFL2eFessd 
os eNO ‘ 
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09/28/76 INPUT LISTING AUTOFLOW CHART SET - LORAN 3 GROUP PHASE FILTER 


FORTRAN MOOULE (LEST»NAMSQ? | 
Caro NO eee CONTENTS eoee | 
1 SUBROUTINE CONF (TAO) | 
2 IMPLEICET REAL SBEAWH 0-2) 
3 DIMENS LON [HC (60)-GROTOTES, 320, ARR EALIGe2) 
4 COMMON IHO,GROTCT,ARR 
5 COMMON/ PR INT/ IPT 
6 AMAK#0.0 
7 00 L 1253262 
8 A=OABS(GAOTOT( LID? 
9 (ALE APaxD GO TOL : 
10 AMAKea 
vt thet 
12 & CONTINUE 
13 AMIN=9999,, 
14 00 2 [#252062 
is A=0A3S50.5-DABSEGROTOTEL fds amaxd) 
16 TF(A.GE AMIND GO TO 2 
a AMINEA | 
18 f2=1 
19 2 CONTINUE { 
20 WALTEC LPT, LOCO) IL : 
2 WRLTEC LPT, LOOLDI2 | 
22 TaOs=SeCL6-f1t - : 
23 RETURN : 
2 1009 FORMAT(SIK,*MAK VOLTAGE PEADING AT Ft, f2) i 
25 LOOL FORMATI// 50K eI VOLTAGE READING AT T'el2d 
2 ENO } 
: 
; 
Program Listing A-42 Frequency Filter Lock-On Confirmation | 
j 
: 
09/28/16 INPUT LISTING AUTOFLOm CHART SET = LIBAN 3 GRUUP PHASE FILTER i 
FORTRAK MCOULE (LEST »NAMSOD | r 
eo NO eeee CONTENTS ease | ; 
1 SURADUTING ALUCKECATSLICT,LPCT TERR? ; 
2 OLMENSEON DAT CE De PF EL 2d, IHL 40) 
3 COMMON I4C 
* COMMON/GE TOAT/ERT SEM UF 
s COMMON/ PRINTS IPT | 
6 [PatPctet . 4 
’ 00 2 f+(o,7r 
8 THOUS DRT RC ESS DOL 
q NaMOOLEMO CIS DOS sod e31 
10 OO bo Jets 
nM IOCT#IOC 92 ; 
2 IFCLER® NEO) GO TU I ; 
1 KeLoett-Ldes 
1s CALL AVAL CCT Ke CAT y JehD 
ie 15 1 CONTINUE 
16 & TECLOCT «EC.16001 GO TN 3 
al ir 2 CONTINUE j 
" 8 3 1PCT#l 
19 RETURN 
20 1000 FCRMAT(I6F 8.4) ' 
21 €40 ; 1 
Program Listing A-43 Frequency Filter Data Blocker | 
; 
f 10/14/16 TNPUT LISTING AUTOFLOW CHART SET = LIRAN 3 GROUP PHASE FILTER : 
FORTRAN MODULE (LEST NaMSOI 
CARO NO eee CONTENTS cece 
1 SuBRIUTIVE TIMECTA) : 
2 OTMENSTON [SUF E1289, THOC SOD 
“a 3 DOUALE PRECISION TA 
4 COMMON IH 
s COMMIN/ CETOAT/ 127, TAUF 
6 COMMON/ PRENT/IPT 
? Th#0.000 
8 OO 1 f#1,? 7 
§ 9 Jeteun 
10 Kele2-1 | 
‘ i CALL BCOTALLE Ss Ke 2 eT BUF. E28 eMD 
4 12 1 T88tM-1HO (2) P86 OPC Tbe TA 
13 CALL BCCT ALLL. Te? e 1BUF 128 64d 
146 T8«TReMe, 000000100 
‘is RETURN 
} 6 END 
¥ Program Listing A-44 Frequency Filter Time Decoder 
% 
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10/1 este INPUT LISTING @UTOFLOW CHART SET - LORAN 3) GROUP PHASE FILTER 
FORTRAN MOOULE (LEST WaMSas 
caro nO eeee COMTENTS eees 

1 SUBROUTINE SETCDAT) 
2 OIMENSTON DATEL? | 
3 OC ft ft, 1600 | 
4 1 OaTtIh=9959, | 
s RETURN | 
s END 


Program Listing A-45 Frequency Filter Array Reset 


09/28/ To INPUT LISTING AUTOFLOw CHART SET - (9RAN 3 GAQUP PHASE FILTER 
FORTRANR MCOULE (LEST NAMSOD 
caad 40 eeee CONTENTS seee 
4 SUBROUTINE RYL EOC TR eOA Teh ND 
2 OIMENS TON CLiaL92) 
3 OUMENSION OATELD, TRUFCL ZA DS THO TSO) 
s DIURLE PRECESLIN BART AL G62 eDSIDTOTES, S20 
s ONUALE PRECISION & 
6 COMMON IMO WGROTCT. 
Es COMM INA GETNAT ZIRT, 
8 CNMYON/ PRINTS LPT 
9 300 
10 tetaurint 
it CALE SHFTUCI ete 32 ed 
i? CALL SHFT ACS  Te- teed ® 
13) MATCLOCT-Liesel 22085 75E- Se THICK 
hd J*0 
is CALL SHFTLO S41 e325) 
16 Cat CHET Say i ei e 4d 
uv DATCIOC Td #54. 2208 8 TSE ~ SO THOEND 
ih) FFE THOE 25 9.GE 64192) GO TO LO 
9 TFEEMOOCTOCT,379.NE.12) GO TOL 
20 THOCZSDREHOC 2501 
. 21 AQRETHOC 2S), de CATELOCT-Ld 
22 ARREIMO MZ 5b ,20eCATEINCTD 
23 c THOU 26) TS SET IN CR8OIN AS WO CYCLE OR PEAK PULSE amP PLOT 
2 b LF LMOOLIDCT 92) NE LPN 26)I) GO 1D 10 
2 ASOATEL OC T-Lpee2eOatdtocrpee2 
an QECTHOL 25 PoeOSuRTOAD 
2r 10 CALL EOITUICCH, CAT etd 
2a RETURN 
29 ENO 


Program Listing A-46 Frequency Filter Measurement Decoder 


09/28/16 INPUT LESTING AUTOFLOW CHART SET - LOXAN 3 GROJP PHASE FILTER 
FORTRAR MODULE (LIST .NAMS QD 
Card NO eese CONTENTS eeee 
1 SUBROUTINE QUTOUTETO OAT) 
2 OLMENSTON Dail) 
3 OIMENSTON GROTNT( 3,32) 
t ‘ OLMENS TON [BUF C1241, fHOL SO) 
s OOURLE PRECISION TMeX¥eKV2ZeKYSOSKRYN 
s OOURLE PRECTSION GROTOT 
’ DOUBLE PRECISION TEMP 
8 COMMON I40,GROTCT 
| cy COMMON / CE TOAT/1°T .18UF 
H to COMMONS PRINT/ TPT 
eI DO 3 ft, 1600 
1 12 TF (OATET) .GT.9000.1 GO TO 1L 
{ rf) JeMOOtT 432) 
is 1FUS6E9-0) ga32 
is GROTOT (Ty JdeGROTOT( Tes) eOArC ts 
ity GROTIT( 2, INeGROTOTIZ. Jd OMAT ELD ee2 
7 GROTOTE Ip JPeGROTOTE Ds sPOR 
1 1 CONTINUE 
19 WRITECZLDTO,COATELD oT el pL6OCd 
i 20 3 RETUIN 
a ENO 


oe 


Program Listing A-47 Frequency Filter Output Tape Writer 
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Fil rvonsutindnenaintn ialthaina’ td iahetalctatlecha diate ok onal 


09/28/16 INPUT LISTING AUTOFLOW CHART SET ~ LORAN 3 GROUP PHASE FILTER 


FORTRAN MOOULE (LEST »NAMSOS 
Caad NO aces CONTENTS eevee 

1 SUBROUTINE EDI TIICCT,OAT, J) 
2 OIMENS TON [B8UF(128),1HO( 40) 
3 OIMENSTCN CAT(L? 
4 COMMON IHD 
s COMMIN/ GE TOAT/12T, [BUF 
6 COMMONS PRINT/ IPT 
? IF (ABSCOATELOCT-199.GT.10.) GO TOL 
8 1F (ABS(CATCLOCT 2).GT.10.9 GO TO § 
q RETURN 
10 L LOCT=10CT-205 
i IHOCSID=LHCC SOD OL 
12 00 2 J=leto 
“a LOCT={0CT 42 
14 OAT( IOC T-1)=9999. 
5 OAT( IOC T) #9999. 
16 2 CONTINUE 
7 RETURN 
18 €NO 


Program Listing A-48 Frequency Filter Editor 


0928/7 we LISTENS AULOFL OW CHART SET - LOQAN 3 G2TIP EuASE FILTES 
FORrRAK MODULE (LIST, SAMS QD 
Caro NO CONTENTS ecoe 
1 SuUaPOITINE STatt 
é OIMENS TON TRUE ELA RDS THOLGOD 
3 OIMENSECN GONTITE S32) 
a OWE PRECISION’ GPITOT, AMAL 
. COMMUN IHC eGRETCT 
6 COM INSGE TOAT/IOT, LPUF 
? COMEONS PAINTS LOT 
a Of b tet. 32 
9 TEEGATE. EP. EI.07 GI TO OE 
lo GROTITEDe PPeGHOTOTELS EI7GROTCT OS ED 
iu GAOTITE ZS TPGRIFITE ALI /GHO TOTES oD j 
12 GPOTOTE 2. LHS MSQKTEGPOTINOTE 2. EI-GQOT ITE LL bO8 Ad 
= iB) 1 CONTINUE 
is TeGOOTITE Ube. S 
6 weLTECror 9070! 
6 WOLTEC TOT, L060 
ur WATTEL LOT ,LOODPECSOTITE Le ED oleh, 320 
1s wetrecter.tovor 
@ APUTECLOT SLOOULEGFOTITG ela elebed2t 
20 WRITECLPT LOOM IMOEITD 
at WETTECTPT LOOwd THOCSaD 
22 w7tTeCT,Locnmd 
2) 1CO4 FOQMATELI A, 10°) 
7 WOETEC TL COODEGPOT ITER ED oh he 2d 
25 WRITECTSLOISHIGRITOTE 2 LD eh ate 22d 
26 3 RETIION 
2? VOOO FIRMATE LK IHF RS 46 /) 
2, LOD] FIAMATE 2K IAFA LSS) 
29 1002 FORMATE*E GRAND TOTALS OVER* othe" PILSES APE AS FOLLOWS TD 
30 1OOS FIAMATES7 LOK, TOTAL NUMAE® OF FEI“CTED PULSES =*e15) 
ot 1006 FURMATES/ LOK TOTAL WIMSF2 OF MISSING OR SHORT RECTROS © 615) 
3 32 1005 FORMATC PFS. 4) 
” 1005 FIRMATO//* MEANS IN VOLTS TIMES LOD 
“ LOOT FOR@MATI//* STANCARO DEVEIATICNS AdiUT THE SEANS*D 
ob T ENO 


Program Listing A-49 Frequency Filter Statistics 
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09/28/76 INPUT LISTING AUTOFLOW CHART SET - LIRAN 3) SROUP PHASE FILTER 


FORTRAN MOCULE (LEST»NAMSQOD 
Caro NO eeee CONTENTS 
1 SUBROUTINE PLOT 
2 OIMENSTOV [HO(460),1BUFC 129) 
3 DIMENSION GROTOTC3 e329, maRK (325 
4 OOUBLE PRECISION GRoTorT 
= COMMOW IHC,GROrtT 
6 COMMON/ GE TOAT/T RT, 1 BUF 
7 COMMON/ PRINT/IPT 
6 WRITECIPT, 1000) 
9 AMK=0 
10 09 1 f#t,32 
ul TFCAMX.LT .OABSEGROTOTILS IPD) AMKsDASS(GROTITOLS TDD 
12 1 CONTINUE 
3 MX=AMKeL 
ind MO=S5/MK 
15 00 10 I*t,55 
le t=55-1 
uv 1€=#107865 7024 
is TFCMOOC Le MDDEQ.O) IC =IL/ HT 90555 3661089470466 
i 00 20 J#1,32 
20 V=DABSCCRCTOTIL SDP 
2. MARKE SD ALCT3IGL A246 
22 Ka5S5-MD ey 
23 1FCT.NE.K) GO TO 20 
24 MARKO SD #~ 9227464908900 342) 0251658240 
25 22 CONTINJE 
26 10 WALTECT Pr. LOOLd IC, MARK 
27 WRITEC TPT, L002) 
24 WPTTECTOT,100%) 
29 ae Tues 
3o LOOD FURMATETL sey LAETT , KD, LAE TS, SDD ry 
iv 2002 FORMATULE AZ LE ,S2EALSIKDD 
32 1002 FORMATE Ox LAC PEt, nde LOCTs® g¥Khd 4 
33 L003 FIRMAT( 20xX6%2.55-6 SECS PER OLVISI ING 4OKet7L56-& SECS PED Olvisto 
36 NSD q 
35 ENO 
Program Listing A-50 Frequency Filter |/O Plotter 
OF/2H/ 16 INPUT LEISTIAG AUTOFLOW CHART SET - LIKAN ¥ G2NU7 PHASE FILTFR® q 
FORTRAN MOTULE TLEST.NAMSOD 
caren NO eeee CONTENTS once 
1 SUBROUTINE SKIPFL 
2 OIMENSION TAHUFEI2Z48D 
3 COMMIN/GE TOAT/E 27, TAUF 
4 COMMONJER INT ICT 
s i CaLt Ge tRec 
6 GD TOCL eho Ve2deERT 
? 2 wairecier,1000) 
8 3 RETURN 
9 10900 FORMATELK,*ND OO CAPO THES FILE CSREPFLOOD 
19 END 
. Program Listing A-51 Frequency Filter File Skip 
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ee ee ee, era eee eo 


AUTOFLOW CHAKT SET - (IRAN 3 GROUP P4ASE FILTER 


09/28/16 INPUT LISTING 


FORTRAN MODULE (LTSTsNAMSOD 


CONTENTS 


Caro NO 


SUBROUTINE PLOTLCAVGPWRD 

OIMENS TON DQI( A192) 

DIMENSION ICH(1L28) 

DIMENSLON MARKC 23 LHOCGOD SE TETEZD 
ONUBLE PRECISION AVGPwa (2), GaP 
DOUBLE PRECISION APRIBL94 42 DeGROTOT ES s 324630256421, SAK, BM IN 
QOURLE PRECISION TOTP Re 

COMMCN IHC, GROTOT, APA e001 SP 

OaTA wARK/* e / 

DATA ITIT/*QQCOe TELITS 

1 COMMON/ PRINT/ IPT 

12 00 6 M=ly2 

13 L=0 

14 TOTPwWR=0.009 

15 OCPWR APRIL MIM e2EADAL2 MIS e2 

16 WRITELIPT »LOOLDITIT&EM De OCP mR 

‘7 NPNT=2ee@THO(25) 

18 LLENPNT /256 


2ertevrune 


a 
- 
° 


19 L2=NPNT -L 102 
20 L3eLbs2 
| 2. OO 2 [a2el2st1 
\ 22 Letel 
23 B(LsM) =. COO 
I 26 00 1 Je delly2 
| 25 Keeley 
| 26 L BEL eMDEACL MD OAIR CK MPO RZOACS (KOL Mee? 
27 TAT PwHeTOTPWR OB (Ly MD) 
, 28 2 CONTINUE 
; 29 BMAK#-9999, 
| 30 BMINED999 « 
| a 0) 3 L#ly256 
| 32 BMAKOOMAXLEACT MD, OMORD 
| 33 BMEN*OMINL CAC TS 4), RMEND = 
, 4 3 CONTINUE 
4 SP0*A4AK-BMIN 
% GAPHL?T/5P0 
” 09 # Leleled 
38 © ICHCT De MARK ELD 
319 00 S [#t,254 
| 49 Ka(BUT MP -AMIN) CGAP OL LOL 
] 41 ICHUK De Maa (2) 
| 42 WOLTECEPT LOGON Ty ICH 
| 4) ICHUM D2 MARKOLD 
| “os 5 CONTINUE 
“5 AVGP we OM) @TCTPMOBLEZNENT 
“6 WITTEC LPT pS LOO2DAMT HMA Ky TI TP ee AVGP WC M) 
“7 6 CONTINUE 
aA RETURN 
ico “9 1000 FORMAT (IK, TS, 12401) 
$0 LOOL FOHMATEOTL Bx, PT TaL NC Pwee IN IK Ale ® H*FF Tow) 
st LOO2 FVAMATE LH SPMIN WELLE PPS LOG yoke PMA VALUE 88 SEL OG OK 
$2 TPT ITAL AC POWER HLF Tee Sky MAYS AC POWER # SF TLS) 
* ENO 


Program Listing A-52 Frequency Filter Envelope Plotter 
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09/24/76 INPUT LISTING AUTOFLCw CHART SET 
FORTRAN ROOLLE {LEST ,NAMS QD 
caro NO CONTENTS 

1 SUBROUTINE HIST 
2 TMOLITET REAL@R(A-HyIT-2) 
3 DIMENSION AFR( 819452) GPOTOTES, 32b-THOCSOd 
bd DIMENSION JIS (20420), 18UCK2 620) 
s COMMON THO GROTCT,aARe 
6 COWMON/ PR INT/ IPT 
7 S$1G623* (CROTOTI 2,11 deGROTOT (24129) 
8 SIG=SIG6/6 
9 IS 1G=5 
10 AINC#2¢S1G66/20 
a BNOO=GZOTOTEL»LLI-SIG6 
12 UBNOQ*GEOTOTEL LL eSIGS 
3 BNOI=GEOTOTI Le 129-S1G6 
4 UBNDI=GROTITE Le L20+S1G6 
15 CALL CLEAR(JIS,S00) 
is CALL CLEAR(I BUC 640) 
uw fa*0 
18 ta+o 
19 raQg*0 
20 avQ=0.0 
au avi+0.0 
22 Tie2e*tO25) 
23 00 20 [lett 
24 KFLG#O 
25 LFeGed 
26 On 3 K*le20 
27 AQHAING #1 K- 2d eH NO 
24 TFCAP20 1, Ld.GE.8Q) GOTT 2 
29 so Ine 
30 2 PECMPAC TL LPL LE. CAQearNC HD GO TOS 
ua 3d CONTINIE 
32 ® TCrtort 
3 KFLGOL 
oro $90 7 Ue1.79 
6 Areaincedc-Lieswe 
v6 fF (e220 1, 27.GF.atd GO TO 6 
u cn To A 
BLJ © LF CATED 2DLE. CAT OAIMC ID GOTO 9 
319 7 COMTINIE 
40 6 tartacl 
ai LFLGEL 
42 LFERFLG OM FLG-LILD 620610 
a3 Tasto-t 
“6 tasta-t 
“5 Tags {Agel 
46 so 19 20 
aT LL STSEK eh de SES EK LD OL 
a3 avGeavues ott td 
ao AVieaAvi oder (ty?) 
so THUCHLe MPR EICO Le KDOL 
st TRUCE? LP eEAUCE2LdeL 
$2 TONT*1ONT 01 
Ss) 20 CONTINUE 
se L2*2e* (HOC 75)-1@14-faad 
35 avaravo/t2 
So avieavi/t2 
ST we TTEC CPT, LIOLD 
SA 59 39 tebeco 
37 x MOTTECEOT LCOOD CILSE ZL bed) esrb 2dde LBUCELS 2E-Ed 
60 wOTTEC tor, too2) 
6. WOLTECTOT LIOVPETAUCE ZK De Reke 200 
62 Me CTECLOT IOSD EN 
6) AP UTECTOT SLOLIISIG 


n 100) 
2 1901 
LB) 1092 
™ 1093 
1s tous 
7 100s 
” 1006 
cL toor 
9 oof 
#0 1009 
La toro 
a2 tos 
”» 


AR LFECT PT LOLODISTG 

wOTTECEPT L205) 10 

MOUITEC TOT, LOOK TA 

wRITECL PT LOOTI IAQ 

we LTEC LPT. LICBIAVS 

wRITECT OT, 1009) AVE 

Qe TURN 

FURMATC LOK 2005. 4K 1 5e/7) 

FOQMATAPL SOK, TLE-T12 HISTOGRAM, 49K, * TLL OSD 
FOAMATE PL Seu, tTL2 DISTRIBUTIONS) 

FORMATC IA Ke 201S///90 
FORPAT( OS K,* TOTAL NO. PATRAS ef 
FTAMATCSIRe*TLE ENLT PAIRS etter 
FORMATIS3A,*TIZ EDIT PAISS ete hed 
EQAMATCSS HST TUL GL? EDLT PALES eh aloes) 
FIAMATOSIK, ITLL AVERAGE eh Teed 
(S3K,9TL2 AVERAGE Freed 
FOUMATEOSIASTENT TING LEVEL Zot SIGMAt esd 
EC QMATCS3KS*TLLeh2 EOTT SIGMA #*eF Tea) 

eNO 


etSesd 


- LORAN DS 


Cane PHASE FILTER 


Program Listing A-53 Frequency Filter Histogram Printer 
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SR Sa a Ree B+ 


09/2 8/ 16 INP UT LISTING 


FORTRAN MODULE 


CaRD NO 


Sane utune 


. 


AUTOFLOW CHART SET 


{LIST »NaMsod 


10 
iB) 


12 


20 


at 
100 


1090 
toot 
1002 
1003 


CONTENTS 


SUBROUTINE TOA 

IMPLICTT REAL*B8(A-CyE~HyO-2) 
DIMENSION DATEB1I92) 

DIMENSION THO(GO) + GROTIT(E Ss 32d APR (BL IG 2), TITLEC 
OIMENSION [BUC(24) LIKE (24) oMS2K(2), 5BUC( 26) 
OaTA TITLE/*BEFORE e,¢ attee *4 
CATA MARK /* ° ‘ 
COMMON IHO,GROTCT,aR2 
COMMING CAQO/ ANG STORSNOES 
COMPON/ PRENT/ IPT 

CALL CLEAR(OAT, 41925 

O09 10) J=le2 

30 TOUR sP)ed 

AMAX2-5595, 

ANIN=9999. 

$100.0 

AvVG=0.0 

TCNTeo 

CALL CLEAR I BUC 24) 
NPNT= 2° e@1rO(25) 

OO 1 L=L,APAT 

TF (mart) 667.9990.) GC TOL 
TCATR(CAT of 

FO TO oT ded 

UF (62201, 2).6f.0.9) GO 10 2 
a°a(tyt apectyt) 

BRACE, 2 De -APR(L 2d 

DATO UPROAR TAN ZESPR OL LEY SO9G F,2EdO1. 591569 
TFCOATOLD CT AMAR) AMAR =OATO(LD 
TF CDATOLDUTSAMEND AMENDDAT OED 
ave ANG EU ATCES 
STI*ShreOariLpes2 

CONTINUE 

AVGRAVG/IONT 

STQ*STOALONT 
STO*OSIRTESTO-AWee2) 
AINCeCAMAK~AM1NI/26 

TCNT #0 

me $ feleSone 

TF crardt) Gt .9990.) GN 19 § 
TIONTSICNT ot 

OD 3 «el, 24 

BUC «aME Re me apur 

TF COATCHD LE.AUCD GC 19 4 
CONT INJE 

Kanel 

PRUCt Cd EAC Cm deOL 

COMTINIE 

MAO 

OF 10 Teles 

TFCERUCCLD GPA) SSeetHuUCOdD 
SD VW Fetes 

SVICCT DST MUCEL DOSY. sterol et 

ME TTECEPT OODPTETLE CD 

0) 29 fel.s0 

OO 12 wl 4 

LINE Up eae 

TRCESL- TP LEQ SMC CRD) LINE Cr pew arr (27 
CONTINJE 

WE TTECT OT, LIOLILINE 

CONTINIE 

MELTECEPT 6002) [8UC 

WITTEC LOT LIC AWGeSTOSICNT 
OF 21 fel ehOnr 

aoarenateny 

TF (0095 (A CAT-avVG).GT. YSTOR OATELI 29999, 
CONTINUE 

CONTINUE 

METAS 

FOGMATU LE S3Ee TOR SeANe*® EXITING (380 CHCLED* S/ 
ECOMATE SC 28C SKA UED 
FOOMAT (SK SIS y 7) 


FJAMATE 26 Ky AVERAGS TON O85 63,4 6% MICBOSECONOS Oke 


LES a4 TNT. OF POINTS 9851S) 
ENO 


~ LOPAN 3 GROUP PHASE FILTER 


$.% 9A of, 


Program Listing A-54 Frequency Filter TOA Distribution Plotter 


nan cate ttimtin: al a el 


= has = 


a ee, ewe eee ae Te 


1 men a 


097 28/ To "INPUT LESTING AUTOFLOw CHART SET - (ORAN 3 GROUP PHASE FILTER 


FORTRAN MOOULE (LIST. NAMS QD ‘ 
Caro NO eoee CONTENTS eeee 
1 SUBROUTINE ame 
2 NMPLICIT REAL CB (A~Ce&~H,0-2) 
3 OIMENSLON [HDU49) -GOOTOT(C3, 32) eo AAR(BL IG e2deTITLECZI-OQIESL92) 
* DIMENS TOW [8UCE 26D-L ENE C26) MASK 2) JBUCE24) 
5 DATA TITLE/* BEFORE "eo AFTER %7 
6 OATA “aRKs* ,ceeeets 
7 COMMON IHC eGACT OT, AFR,O0L 
8 COMMON/ PRIAT/[ PT 
9 09 100 Jele2 
10 GO TOL8eFded 
ut B AMAX=-9999, 
12 AMIN=7999., 
Le) 9 $TD=0.0 
iad AvG=0.0 
1s ICANT #0 
16 CALL CLEAR(T BIC .246) 
‘7 MPNT s2e ef HOE 25) 
16 OO Lb T=beAPNT 
19 lf (COQ1(1).GT.9000.9 SC TO 1 
20 LONTeIONtelL 
24 GO TAHT Hes 
22 & LTFCOQLELD -OF AMAR D Awake DOTS 
2 CECOITE TD LT AM INDAM IN ROUTE ED 
24 T AvG#OILCLdeAvG 
25 SMOedILdLaredestn 
26 1 CONTINUE = 
27 AVGSAVG/ECNE 
2A STO<STISICNT 
29 STO#OSQeT ISTC-avGee2) 
30 ALNC#(APAR-AMINIS 26 
34 {CNTed 
32 00 5S Tele SPNT 
33 te COOrtt).6f.9000.9 Gn ro 5 
oy LCNTOICNT OL i: 
35 NO 3 Ket, 24 
we BUC east te tw 
uv TECMOTE TD LEMICE GO TD © 
33 3 CONTINUE 
39 Rakel 
40 & Pau he del BC Cepet 
“1 S CONTINUE 
“2 A460 
a) Of 19 Lebede 
as LO TRULAICCLD.GT. AAD PAxe BUC OE? 
“6 OO LL Leber 
“6 TL PHOT ELD RP PICE1D OSD sa OL. 
a? WOTTECTON LOOP TE TLE CSD etn d2od 
“A 09 20 teleS0 
ay On t2 mel ae 
so LEME CRS eA wT 
st. TECESL~ EP EQ. COD CINE Cm Des am dQ) 
$2 12 Continue 
ss WRETCC TOTS LOOLILENE 
“ 22 CINTINSE 
“~ ss wOTTECter loose fHuc 

36 WHITECTOT  LIDSPAVG eS STOSICNE 
sr DO 21 Eel ehPnr 
34 Qreogretr 
sd CECOASS EDT H-4UGT GT. STOR DOLE 9299, 
60 21 Continue 
er 100 CONTINUE 
62 TECIMOC 2 DLE QE LD CETION 
e) GOEOTOTCEIMOL2 ADD ea, 
os GAITITE 2, EHEC 26d eS TD 
6s aerues 
on LOO FHM STE ID SSE, AMPL IMIDE % Abe EDITING oe 2Hy MCSAMPLE "ede * dt eZ) 
or HOOL ED@MATOC SE. 24 aun : 

: 6 1002 Foewarise,26i5,/3 

> Cd DODDS FQAMATE QTE AVEPAGE AMOLETUME © FAS e® VUTS Om, Se. AMPLITUDE 
nD Let eFA eRe ND. OF PIINTS | 8ES 67) 
mn eno 
Program Listing A-55 Frequency Filter Amplitude Distribution Plotter 
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oeszaste taeut CESTING AUTOFLOW CHART SET - (GRAN 3 GROUP PHASE FILTER 
FORTRAR MODULE (UIST. NAMSQD 
Caro NO eeee CONTENTS eeee 
1 SUBROUTINE EDITE(AVGP mR? 
2 IMPLITIT SF aceata E-H,0-2) 
3 DIMENSION 2625652) -AVGP m2I2) 02116192) 
* OUMENSEON 252081942 eEHOLGODSGPOT ITE 3, 32) 
s COMMON IHC. AOTCT APR eOCI OB 
6 COMMON/ PE INTs IPT 
7 NONT#2ee[r 7425) 
8 WRETECTPT «to 90) 
9 00 100 f#1.2 
: 10 awatsd 
iu LL=NPNT /25% 
12 09 10 J=be2%4 
3 TECBCSe ED LT EMOC2TISAVGPMR ITD) G) TO 10 
ind L2"(J-3deL led | 
5s L3eL2e5% 1-2 
16 TFCL20T. 3D L283 
uw DD S Ket2-t3 
18 PFCK.GT Atte) 67 19 S | 
19 (FCARICK,ID.£ 2.0) GC TOS | 
20 aerating fbeC.2> 
at TONTHICAT at | 
22 S CONTINUE 
23 12 CONTINUE 
24 Tl eloet | 
25 WIITEC LOT LIDIDICNToEE | 
26 109 CONTINIE 
2 RETURN 
24 1099 FISMATE tbe SSee TALL JUTPUT DATA SGA te tEML OMEN THES SECTL ING Se | 
29 WOO Ke AL LL FF FPEQUENCV eS eSOKe FILEFOFOt ASIA 
q 30 LOCL FOFMATC SON, .* POLATS FELT EF ZED SU OW Tee leest 

x to 


Program Listing A-56 Frequency Filter Notch 


| 
OWs2as IA teeut COSTING BITOEL OW CHAST SET = LD4AN 9 GAO UP Pease ECILTES | 
] 
FORTBAR MOCULE CULSTNAMS OD ] 
| 
cae NO eeee CONTERTS eres | 
1 SUA IITING wysee 
2 PPPLICEY SF a, edt a-C f= OW- 0d 
, TUMERS LOG Ge EGG pCR TECT RSS S2b BAG EHE DG ee be DIEUALIDD 
4 COMMON PO, DE OTOT, a4 5098 
s COMMONS PRINT SLOT 
6 NPNTa2® erm 
7 ov bt tet 
8 Ouldtd. 
9 1 CONTINUE 
1o RETURN 
uw END 
~ 
Program Listing A-57 Frequency Filter Data Selector 
097247 16 INPUT LESTING AUTOFLIM CHAST SFT = UTRAN § GOI? MHASE FILTER 
FoarRas MOOULF (LEST NAS) 
cay KO eee CONTENTS eeee 
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Program Listing A-59 Group Phase Library Raw Data Tape Reader 


macko 
Saves 


AUTOFLOW CHaae SEL ~ LIRSN Ganue PHASE LIBRARY 


CONTENTS 


£4, GA, £0 


EXTENDED SAVE MACRO 
THIS MACRO PROVIDES THE FOLLOWING? 


ISSUES @ STANDARD SAVE MACRO (CONTROC PROGRAM SERVICES) 
GENERATES AN LA-woRO SAVE ARFA 

PEPRORMS LINK ACE JF HIGHER AND LIwER SAVE AREAS 

INSE2TS CURRENT SAVE AREA ADDRESS IN REG U3 

DOES NIT ALTER GInee PEGISTERS 


WHEN USING SAVEM, EXIT May BE ACCOMPLISHED BY? 


L L3estt dd 

RETUGN CNYME 

cere Exe 

SETC *KQKQ*, PESYSep ae 

cs oc NEED & NAME TQ COUNT FROM 

SAVE FABRA JUST AS IN C.P 65. MAUUSL 

ir thel3 HOLO 13 (HIGHER SAVE AREAD 

CNOP 044 ALIGN SAVE AaEs 

BAL LD eB-EKOS*TOCLSE JUMP PAST SAVE AREA 

oc LaF eo SAVE AREA. cel 3 POENTS MELE NOW 

ST 13,8011) PUT LIwER ADUZESS IN MIGHE? SAVE Ay 
st thew Ct3E AND HIGHER ADDRESS I LOWE? SAVE Ay 
t LivOst ta? RESTORE FROM HIGHER SAVE AREA 

MENO 

SPACE 2 

USAGE: 


CALL GETREC, OR 
CALL GETLe 


WHEZE USE® HAS DFEINED A LABELLED CCMMIN GENO ATS ERA DATA~ 

RHERE E22 (S A FULL WORD AND DATA HAS & LENGTH JF 

AT LEAST "LRECL* BYTES. 

Space 2 

csecr 

equ ° 

ENT@Y GETL 

space 2 

SOVEM E14, 1L2Iee® 

ue 10915 

USING GETREC »10 

L LL,=v(GETDaT) 

USING GETCAT ALL 

ta 12,0C8 

USING IHACC3,12 

space 

Ab OCBOFLGS+x* 10" 

80 Ger JF OCA OPENe a0 RECIFO 

space 

TREN OCE, INPUTS 

th CCHOFLGS eK" 10° 

8 nano SET ERROR FLAG TO NC “ORE TAPES 

space 2 

EQ . 

Lh ERWKS) EXLG*LY INITEALEZE FTO GOOD 1/0 

GFT CCB, INCY 

SPACE 2 

flu ° 

u V3.4) 

RETURN [heel 2det 

eect 

EQu . 

st O,ADFSSI CONTAINS ADDRESS OF STANDACO 
STATUS TNOICATIAS AND & OLSPLACEMENT 
VALUE TO REACH THE CCw 

st t.ageuca CONTAINS AGO2ESS OF CCS 

wc Caw) oeerar2* SET TD 179 ERROX 

wR 0 Pur 108 ANO@ESS INTO RL 

wc Savi99t20r,oth) IVE INTO LOCAL AgEA FOR 

BEFERENCE 1F CALLING PRIGRAM TEAMINGTES 

wiTH & OUMP 

pe 14 BETURN TO SYSTEM, WHICH WILL 
RETURN TO STATMENT AFTF2 GET 

space 2 

eq e 

wve ERw(4) oe XLG*3* SET TO END OF FILE 

cLOse (oce? 

space 2 

CHANGE DONAME BY [NCREMENTING *TNPUTL® TJ *tnpuT2® ETc. 

uo %, HO" 

tc 4 OC 300NEMO% 

an GerHele 

cH 4yaHt 250" 

ac SOLESS 

ic 4,00 80 NAM +S 

La Seanett 

STZ 4 yDCADONANES 

1c 4, <K°FO" 

src 4, OC 30 0NAM OS 

RIN 

space 2 

€Qu s 

uve Eawieresncees* SEF TO NI OO CARO 

a RIN 

Space 2 

Entay TO CLOSE A FILE ANC INCREMENT THE JONAME, 


SETS UP THE FASE VEGISTER SO THar IT WILL HAVE THE 


Saves 
ENtey 
USING 
t 
pene 
USING 
L 


TMT TAB ASR STIR ET RRA TER BIE + toe « 


UE THROUGHOUT THE CSECT. IT THEN PERF DONS 
PRICEOURE AS THCUGH AN ENO-OF-FILE wWEXE ENCOUNTE®ED 


Clee l2bee® 
CFILE 

CFILE StS 

10, 8( OE TREC? 
is 

GETPEC FLO 
ALesvGGercart 
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Ger 
cer 
cer 
Get 
oer 
cet 
SET 
cer 
cet 
Ger 


cer 
oer 
Ger 
Ger 


0009 
9005 
Q0l10 
9015 
0020 
9025 
0030 
0035 
0040 
0045 
0059 
0054 
0059 
0065 
907s 
0080 
0085 
0090 
009s 
2199 
o1os 
old 
Outs 
0120 
0125 
01390 
O15 
0169 
oes 
ose 
O1ss 
0160 
0165 
oro 
OL7Ts 
oko 
0185 
argo 
O195 
0290 
0205 
o2t0 
o21s 
0220 
0225 
0230 
0235 
9240 
0245 
9250 
0255 
0260 
0265 
o270 
0275 
o2a80 


o281 
0282 
0283 
a2ae 
0285 
0286 
0287 
0268 
0289 
0300 
9305 
0310 
o315 
0320 


o32t 
0322 
0323 
os2e 


0328 
0329 
0330 
1330 
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110 
tu 
112 
113 
iis 
115 
116 
Li? 
118 
19 
120 
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122 
123 
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INPUT LISTING 


AOFSSI 
AOFOCB 
SAVIOB 


2cA8 
GEToaT 


ERw 
TNCY 


La 

8 
EvecrT 
oc 

cc 

oc 
SPACE 
oca 


SPACE 
OSECT 
os 

os 
0723) 
ENO 


@UTOFLOW THAPT SET ~ LORAN 


CONTENTS 


12,0€8 
£0F > 


Fror 

Ftor 

2cx*00° 

2 
OSCRG=PS,MACRFE=(GMD,OONAME=INPUTOLY 
EQOAD=EOF » SYNAD#ER ROR -EROPT =ACC 

2 


F 
DSORG=(QS) -DEVOaITAD 
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1332 
O33t 


0333 
0335 
0340 
0345 
0350 
0355 
0360 
0365 
0370 
0375 
0380 
0385 


09/28/16 INPUT LISTING AUTORL Om CHART SET - LORAN ZROUP PHASE CIBRARY 


ASSEMALY MOOULE {LEST»NAMSQ) 
cago NO sees CONTENTS eeee 

Lt @ FORTRAN CALLABLE PCUTIAE TO COACATENATE Tao FULL WORDS,» SHIFT 
2 © LIGICALLY LEFT OP RIGHT 4 SPECIFIED NUMBER OF PLACES, AND STORE 
3 ® {O-ORDER AND KI-ORDER RESULT BACK IN ORIGINAL WIROS. 
* @ CALLING SEQUENCE: CALL SHFTL Oe A,C,2C) WHER= A AND B ARE OR EGINAL 
5 © TwO WO20S, C IS NUMBER OF SHIFTS &ND SIGN BIT e FOR LEFT SHIFT 
6 ® AND - FOR RIGHT SHIFT, AND RC IS RETIAN CODE. RC=l6 IF ICL 26s. 
’ . 
8 . 
9 SHFTL CSECT 
10 SAVE Clo, t2dee® 
aL USING SHFTL» 12 
12 tr 12.05 
13 st 13,SAVEe4 
14 ce 2013 
15 ta L3 SAVE 
6 st 13-8(2) 
7 tr 201 paraMETER LIST POINTER 
18 La 10,0 SET R8C2C 
19 . 
20 © PICK UP FIRST AND SECOND PARAMETERS IN REGS 6 AND T 
rat t 6,010.2) ADDR OF A 
22 t 4,006) o) 
23 t 7402) 
24 t 5,007) 8 
25 . 
26 hd PICK UP 3RD PARAMETER ANO GET NB OF SHIFTS 

q 27 t 8,802) 
28 t 9008) Cc 
29 Lz 9.9 POS vatuEe OF C IN 29 
30 ure 9,9 SEE tf AC. SHIFTS = 9 
3L BE Exit 
32 c 9y=F 64" SEE IF INS. SHIFTSI > 64 
33 BH GR64 ves 

z 34 ™ 018) .x*80" tS C NEG 

‘ 35 BND POS 
36 tR SROL Sel SHIFT RT SINCE C IS NEG 
a7 ‘ act SLR 
34 8 store 
39 Pos SLOL 4el SHIFT LEFT SINCE C POS 
490 Bcr 9,POS LOOP FOR ALL 
4. . 
42 @® STORE RC ANDO EXIT 
43 STORE st 4,006) A 
44 st 5.0(7) 8 
“5 ExIT t Ite 1L2e2) aOOR OF RC 
46 st 10,0011) STOCE RETURN CODE IN 20 
47 t 13,SAVEeS 
48 RETURN (14,12) 
4? . 
50 SROS ta 10,16 
SL 8 exit 
52 Ge 

~ $3 SAVE os Ler 

S46 ENO SHFTL 
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ASSEMBLY MODULE 


caro ND 


- 
CLanrrrune 


IN@ UT LISTING 


(LIST, NAMSQD 


AUTOF LOW CHART SET - LORAN GROUP PHASE ( LBaaRY 


CONTENTS eee 


FORTRAN CALLASLE ROUTINE TD CONCATENATS Tr) FULL WORDS, SHIFT 
ARITHMETICALLY LEFT OR RIGHT A SPECIFIED “UNSER OF PLACES, AND STORE 
LO-ORDER AND HI-ORMER @ESULT BACK IN ORIGINAL WORDS. 

CALLING SEQUENCE: CALL SHFTA( AS ¢CeRC) WHERE A AND B ARE DRIGINAL 


TS NUMBER CF SHIFTS AND SIGN AIT ¢ FOR LEFT SHIFT 


AND - FIX RIGHT SHIFT, AND RC IS RETUIN CODE} RCHL6 TF ICI >6%. 


. 
. 
. 
. 
© Two wOROS. C 
- 
. 
. 
s 


wETA CSECT 
SAVE (14 sl2dey® 
USING SHE TA, 12 
tr 12,15 
st L3SAVESS 
Rr 213 
ta L3,S4VE 
st 13,802) 
ur 2st PRAAMETER LIST POINTER 
a 10,0 SET RCO 
. 
© PICK UP FIRST AND SECOND PARAMETERS IN REGS 6 AND 7 
t 6001062) a00R OF A 
L 4-06) 4 
L T4102) 
t 5,0(T) 8 
. 
© PICK UP 320 PAQAMETER AND GET NB OF SHIFTS 
L 8.812) = 
i 9,008) c 
tee 9,9 POS VALUE OF C IN 29 
itr 909 SEE IF NO. SHIFTS = 0 
BE Ext 
c QesF*5o* SEE IF INS. SHIFTS] > 64 
BH GRO YES 
™ O(81,x*80" IS C NEG 
BNO Pos 
Le SRO 41 SHIFT RT SINCE C IS NEG 
act FAR 
8 STORE 
POS SLOA 41 SHIFT LEFT SINCE C POS 
BCT 9,905 LOGP FSA ALL 
. 
© STORE RC AND EXIT 
STORE st 41016) A 
ST 5.007) 8 
ExIT L Ale1202) = ADDR OF RC 
st 10,0011) STORE RETURN CODE IN RC 
t 13,SAVEes 
REFION (14,12) 
. 
Gees LA 10, 1é 
8 exit 
. 
Save os 1aF 
END SHFTA 
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09/28/16 INPUT LISTIAG auTO*LOw Crear Sf - (Fan Gague PHASE Lf sasey 
ASSESBLY MODULE UIST. Namt5ay 
ca 4a weer 
1 DFETA t r OFT 1 
2 © pFeta ne Ober we * ’ ISTON fast & WISP TePNSeoas, THs ofta 2 
’ . VERSION 15 elt foeteay iv. OF ta 3 
. . us CALL DEF TACA. =, 008T) oth 
$ e aw sPLE PF ECL FLOAT IN POINT IN@UT as OQut®ul Seear, DFTA 5 
* ® 2 1S Ofmenstow gees FOR “PLEX CALCUL ar {OnS. of ta o 
’ . IN THIS CASE QOL) TS FIRST REAL PART, 812) TS FLaSt ofra of 
4 ° {$aGINARY GAPT, 203) 1S SECOND REAL PET, O€e) TS SECOND OF TA B 
% * teaciuaey eset FIC. uP (4-1) COMSLER NUMTERS. 9 
to « A 1S DIMENSTON=2 Zoomer AL CALCULATIONS 3 
it ® IN THIS CaSe Sfat values UL INTO BC bbeeces t 
: ® at reedd, SPFICTENTS Can/2ds2s 
is . BN/2=0 MUST BE PUT INTO RESPECTIVELY. 
ed e - tS ©ta PCTAT TATE Of VALUES 
‘g " . IN A It 6S us 
- 6 . xoer & EIKED pF 
? - 1F MOP T*O, [NAUT 1S tre 
ts > {F aMeTs2, [NeuUT IS C ts 
1 . IF weheetzepreysus, FISD © COMPLER DATA 
, 20 . COMPLE CE FQLaLES Te 1S THe 
2t . Foarsa taay 
q 22 Ly FJD = SUMEC CR Domest sendts FeO de 
a) ° BHEGE THE SuHMaATICN 1S so)! eeNek 
as be ctmd © Lenesuvee i) et te 
25 . wrERE THE SueMaT -1 
; as . TF KOPT=3, pRPUT ES 
: 2? . FCURLES TeeNsrce 
{ ; 2A - IF KOPTO3, INPLT TS 
q 24 . IS IvERSE COMPLER Enus Tes 
30 * TE MCd_ ND COMPUTATION wit TAKE 9. ace. 
b 31 #615 USE BRANCH JEGISTE? AS BAS 
32 ase 
33 Se3 
4 as 8 OFerat 
35 oc RLLTO5° 
ie oc c*orrrae 
% 5 7 BEFTAY A¥= (14,12) SAVE BEGISTFESS IN SAVE APES 
~ 2 38 pane 15 
‘ 3 39 BAL? 100 
q a USING ©,!10 STANCARD PP 1G7S4 BASE REGISTES 
,, at te 3413 OLD S4VE aFEA ADURESS 
42 ta 13,8988 N&a SAVE APES ANDPESS 
a 43 st 13.810,3) ADOZESS OF % SAVE A965 (NTT OLD BREA 
4 “4 st 344C0, 13) ADDSESS TF SLO APEA TO WCED 2 CF NEW Bees 46 
4 45 MVC | AMKL7D,OCLD «= PARAMETER ADDZESSES 45 
7 46 “vr eo 
L oT L aNOR ESS OF 4 TT a6 s? 
t 48 t - “8 
49 € “6 
so ac 1? DFETF20 Exit 2© "<3 50 
st DFFTA2 t 2-40 apogess JF A St 
$2 5 202F ADOZESS MF 4 ane 62 
2 t Sates s3 
s+ eve ener ,oter STORE «ner 56 
3s “ve DEFTAL AOS FILL TNSTRUCTION €O2 © [vas Naay PSeT 55 
{ i ua Set 1 ro a5 56 
$s? SLX Sy 4KaE (Qeeupes FOX B BYTE @0805 TO BS s? 
3 b +4 st STOZE NUMSE2 ELEMENTS IN A SPRAY 58 
% sy ve ss 
°? sea N/2 FLOATING POINT FOR FT FACTOR 690 
P s st iabel ot 
; 5 wa AM FCO 8 4¥Te wR 62 
3 st 83 
as $o° 64 
— “as rs ‘5 
bid s? FLOSTING POINT N/2 6é 
q P er Le eT 
oa saa 68 
e% s NO FOR @ 4yTe W805 69 
ro sr TO 
" c 1 
a el t CHR2ENT ANDZ=SS OF 5 TO AB 72 
ry * DP NCH GENEIATE § STURAGE TF EQuse ™ 
y ™ ce 7% 
s rs ac WO STOSAGE TO PELEAS® IF OLON=9. TS 
bed t Te 
Lad ‘ ? 
i) a COFRECT S ADNZESS TD RELEASE 18 
" FERS MM EN Beet Od. ae tlt Lad 
wed oreTay s’ s.7L08 New “> 80 
a is o.5 al 
r $24 Ge? Re 
. ry O. Fear N/eeol ~ B BYTE wrers - BEQUIZFO FOF § ay 
st O.ckNS as 
GEtwath Fei vel ol as 
, a Jet a6 
s doce tae S aOceess TO a) ar 
st 2.588 ae 
ua 8.7 a9 
s e,se tar L-2 £08 4 BYTE woROS 90 
WwW uo Op sO SL SESG2HSESSITT IIE OD 
92 090 CoEN 
. 93 sto Oe! 2e>n yw 1 at 
bbd ia 1.486 anGueenr 1S at 
$5 t US.¥OSIN as 
* Race 16.1 Sinceri 19 OS 9 
‘ s? STO 0-98 er 
: nt 9a ta Laas ARGUMENT TS 8t 98 
a9 t 15.9008 a9 
too Bare 14.15 castath to oC 109 
tor sto 0,0c iol 
102 soe 0.9 0 fo FRO 102 
‘ tos $f) OeSe* oro sap 103 
§ Los $0846 104 
105 Lo O.*0°L.Eeo* 105 
te 106 una 6.0 104 
4 tor STO) 0.504 1.80 TO StN2=1) 107 
ton io O.0S toa 
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Cako NO 


lov 
iio 
uta 
iz 
au 
Lis 
fis 
te 
i? 
iis 
ii 
120 
tat 
122 
123 
126 
125 
126 
t27 
126 
129 
130 
rat 
132 
133 
lh 
135 
136 
a7 
134 
139 
149 
Lat 
42 
Las 
146 
145 
146 
is? 
148 
149 
150 
1st 
162 
153 
156 
155 
156 
is? 
156 
159 
160 
1é6t 
162 
es 
144 
1465 
1A6 
ler 
loa 
169 
170 
um 
172 
173 
lls 
1s 
116 
177 
1738 
179 
1a0 
tet 
a2 
183 
Las 
145 
186 
187 
1848 
139 
190 
iat 
192 
193 
196 
195 
1% 
tar 
198 
199 
700 
201 
. 202 
203 
7% 
205 
206 
2cr 
208 
209 
210 
zu 
212 
au 
$ 216 
21s 

26 

217 

zie 

219 


as 


Ba a es a 


INPUT LESTING 


OFFT AS 


OFF TAS 


OFer ae 


JFFTA? 


oFeTaR 


OFFTAG 
DFFTALO 


OFFTALL 


DFE TAL? 
OFFTALS 


OFFTALS 


DFFTALS 


Ome 


OvSeld 
S424 

14,8 

1545 

1563 
1S,8Ftae 
O.8 

Bete 

0-0C 

OF 

20 

a“ 

S06 

Oy? 

610 

0.05 
26Sta) 
Osteo) 

9. 14, DFFT AG 
KOPT 3 eX*02° 
e.OFFTALT 
2eOFFTAT 
es ee 
SEN 
CFFT BL 3 eS 
5.3 

14,l6 
L5eN 

0.8 
O,Ata) 

2 aca 
C4 

24 
CoAta 
2080809) 
9614 DEFT 46 
1a 

1? 

alo 

Son 

08 

Tel 
Oat 
2,4ea(Td 
SACL) 

be Aeatt) 
12.0FFTA9 
Ovatit 
2.40801) 
eet 

by AeBeTD 
BS 

det 

TA 
2.0FFTALO 
1.8 

TA 
166,0FF TAA 
B16 

98 

beNH 

146 

Tet 
Ty=F4Bt 
TeMH2PL 
150N 
Os 

1669 

Leek 
TeNHZ2PL 
14,15 
StF tBe 
Lee 

2,OFF TALE 
Ostet 

” 

Illes 
ILeNo 
2eStity 
Br B DFET ALG 
7 

4.0 
MEET ALG 
202 


Oe 
ZeAestitd 
eeatit) 
26 


eae Re sd 

4 Aust 
eeneaes) 
Ovacst 

2 heats) 
SALUD 

er eemeiy 
liyts 

Se LA DFE TALS 


oe wks + eee ae TD 


AuTOFLOw CHaet SET ~ LORaN 


CONTENTS 


SINCett TO se2d 
INDEX 
INCREMENT 


COMP ARAND N/BOL FOR 6 BYTE WIAs 
ALIGN LOOP IN DQUBLF-woPo AQUNDAIY 
L-f TO INDEX REGISTES ~ 4 d¥TE WORDS 


Fe2 TO Stt-t) 
Feo TO Stl) 


SINC 2°P LANG Tbe S=Oebewee NS EN SUL) peee 


REAL SYNIMES!S LF KO? 
Pur apeay tw BIT REWEAS 


L/N TO FES 
FILE INSTRUCTION F22 - IMAGINSRY PaeT 
b - @ BYTE HOADS — TD ENDER FIGISTER 
2 - @AYTE wIRDS - TD INCBEYEN 
W WORDS Foe COMPRO SND BECIST 

ALIGN LOOP ON DOUSLE-«ceD ATU 


ser 


MULTIPLY A Ae ay BY LPN 
10 J - 8 SYTE wO225 
LTD INTER L - @ BYT= =9RDS 
2 TO INCREWENT REGISTER - 
ws TO eS 

ALIGR LON IN ODUELE-w 240 BTUNDAPY 
compare y To t 

erat ad) 

IMaG 200) 

wear act 

IMAG ACT) - SLL TO 2ESESTERS 

Jets DD NOT Swap 


Swap ATE BND ALJ) ~ BOTH 
REAL AND [rAGINGeY Parts 
N (DR 2eN4) TO «K 

Ke2 TK 

J-e m5 

TF Joe, GD FACw 


J-Kowee 2 Jon TOS 

GO BacK [TF IC=N 

2 - 8 BYTE wIROS - TI FA 
2ro! 

NH TO ID 


NH/ZeL TO 27 - B BYT= w>RDS 


NTO ALS 
ws fo Is 


2et m 1st 

NH/2eL TO (8 

Compase (st TON 

1 10 IND - 3 BYTE «2905 
tro s 

FINISHED TF LSTON 
SUINCH TO FO 


LOR O60 32 UNDP O69 - PUSSTPLE SIGN CHANGE 


IND TO BLL 


SUNQ-INOd TD FRA2 
BRANCH TE IND e1DC#12 

oro 1 

INO=15 (0@ %ef01 TO IND 

CO NCT CHANGE SIGN F25T 190° 
CHANGE SI3% OF SENQ-ENOI 

STORE IMAGINARY PAIT TE CUSEENT w 
STO@E PEM PABT OF CyaRENT & 

310 « 


mer fot 


coeatl) to €RO 
Stwacteth To Fes 

CIAL I-STPACL OLE 19 FBO OF aa 
Coeaiiety 1) Fae 

Steaier To Fas 
CO*atLetiesrea(cr Fa FR2 OF aat 
Aa? TO £26 


€a0 
ateetieaat 19 £22 
Atm j-802 TD 06 
a(molp-sat TO Fe5 
F290 10 
FR2 tO 
£a4 10 
Fea TO Atvell 

Lerst ce veer eot T7 t 
STOP TF molsTon 


165 = 


© Fa IF KOPTSS 


aegistea 


$ 
a eyYTe wIeD 


219 
2ut 
212 
23 
216 
ais 
2te 
217 
21a 
219 
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09/28/16 INPUT CISTING AUTOFLOW CHaat SET ~ LOKAN GHQuP omeSh LISaes 


Caro AO CONTENTS soae | 
7 220 ARLE byA.MEFTALZ «sfor LF sot 
z2t ua Gets (st 1 1 
222 SRA Ont 19/2 To US 
223 saa él tov2 TO 10 
: 228 ac AS, OFF TALL | 
q 256 FF IAIS OTS wer X*O1* fF LCw NADER BIT IS ty ecor = LOK Y 
226 ac 1,0FFTA20 Ext’ (© tow O@deR AIT IS bt 
Z 227 OFFTALY La Sit 2- 6 BYTE wOROS - TOT 
228 ua 168 1 - 6 BYTE wORDS - TO et 
229 t TeNO | 
230 se 78 no-t ro ar 
2 ta S124 3 - & BYTE woe0S - TO INITEAL K 
232 we Bee 2 — 8 BITS woROS : 
233 ry eWN 
4 236 sz 8.9 N-Ke2 TO ENETIAL J 
235 ie lars 2 - GB RYTE wOBOS - TD INCREYENT REGISTER 
236 t LS aNet 
aur sa Set NH-1 TO COMP ABAND SEGISTER 
238 OFFTala LO O,aca) AGE TS FRG 
239 pm) 2eAe BIS) ACkeid FD FR2 | 
260 to aoa ad ACJeli TD £24 ] 
24h to ACJ) TO FRA 
262 49 accbestui 17 B48 
243 so 2,A088) AcKkeli-atsety TO ast 
aes 40 HAORE DD ACJetier(Kely TO ane 
245 s0 eearod acsi-a¢ed WO 28F 
245 STO LBD STGEE £42 
247 $19 20h a0 STOSE Aa 
{ 248 AB F|O 
269 agl TO FR2 
250 agzecn 
251 aateco 
232 ageest f 
253 aptest | 
246 fet To ¢ 
255 NQ-1 TO RT 
284 aamecheakes, TO aw 
2sT aeResi-ableco TO ART 
258 ~aw TO £R2 
259 ast 1h FR6 
269 Aste be 
2oL a0 2eA08) aes ~ hw 
262 so 4 Aeatad -bAleaat 
263 ao Ae Add aatecsl 
26% #0? 0,0 
265 wor 202 
266 HO? ays 
26Tr mOR O64 
268 STO Oat eSeCeseaay TO ated 
269 STO 902,489 o52CRfo-24) TT ACD 
270 $to 4, Aeata) eSUU ZA Leas Ty TO Aeon 
2m STO AP BERD eSeCAALoeaty TO atsery 
272 ~ Rete gt wes 
27 BEL= Dele eNFFT ALA STOP IF K>INHHLD 
27% uo Oraed act) TO aw 
275 pe 20 
275 4 16N Rel TO {NOFK SEGISTES 
2 CLL KT e3.KPO2" 
278 Lis BRANCH LE ROPT ae? 
279 ao Kot 2 
280 sto Awehivel) TO ACL) 
2aL $0 
282 sto Aw-StNetb 7) a0) 
1268. “vi CHANGE DPTIny TO 3 
284 8: Go TC RIT BEvEeSat SECTION OFTE 284 
zes OFFTAL® 40 DFTA 285 | 
~ 286 so Fre 286 | 
2ar Hor oFra 287 
288 “02 o= Ts 288 
249 soe eta 289 
290 sto 25eCAwes( 2) TD Atl oFta 290 
2a sto e580 be- A020) TD ALNOLD DFTA 291 
292 sto o TO at2r deta 292 
293 sto oT a(Ne2) oFra 293 
2% OFFTa20 L CLO SAVE aPea ADDRESS a=Ta 294 
295 us Le t2vl2C19) RESTORE REGISTER CONTENTS OFTA 295 
5 296 MVE 12013) .KTFE® | 
27 ace L5yte Re TURN 
298 a osecr 
299 s osecr 
300 DFFTA csecr 
you FN os ° 
302 at os o 
303 os os ° 
304 oc cs c 
305 st os 0 
306 co os o 
yor FL cc K*6E00000 100090000" 
308 a9 05 3 
309 OLoM oc Froe 
y10 SLs ps fF 
ytt . os 
312 “Mw wht . 
uy no os f 
sis WH2PL cs F 
5 saa es a 
vis «oor 5 fe 
yar AER os ae 
ya 49g oc a*80° 
9 sue o2 acueate 
320 es tor o.0 INST@UC TION FOR © IMAGINSEH PART 
32h ms iNO% 0,0 INSTRUCTION FOR = [44GIN 
322 vosiw oc WtOStNt 
32) wees oc wtocost 
sas en 


ae 
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” ST ae Ae, ee Te oe ak LE 


O9/28/ To 


Foar 


N SOOULE 


caeo NC 


INPUT LEST ING AUTOFLOW CHART SET = LO#AN GROUP PHASE LIBRARY 


(LEST NAMSQ) 


oer CONTENTS eter 


SUBROUTINE PRATIL(L RO. TPT) 
DOURLE PRECISION Al ty 8) 
OIMENSICN 1HD(S0) 

WRETE( TON .2000) 

WRITEC TOT» 1006) 

WRIETEC LOT LOOT) 

WPLTECEPT ¢LOOK) 

WOITE( Ter, 1009) 

WAETE (TPT 5.210) 

wWeTTECTeT, LOLs 

WPITECT ORs 2000) 

wALTECL ON LOL 

MOITEC LPT, 1OL3d 

WP LTEC LOT Ola 

WRUETEC LPT LOLS) 

we rec ter, tore) 

weTTEd@er tour 

WOLTECLET 2000) 

WRITEC TOT sLOLAD 

MATTECT PD LOLG) 

MO LTEC (OT, LO298 

MRTTEC LOT LO21) 

WELTEC TOT 10273 

MPITEC TOP» 10239 

WRETEC LPT 2000) 

CALL DATECIMDEL De LOeEMD 

THOC22 #10 

THOC2 35 <1 

CALL DOTOBIA SIS ID fHOU QI) 

O00 3 K=156 

MRITECTPT 2002V (ACH ID debe Ad 

MRITEC TPT» 2001) 

MeLTEC TOT L900) 

wRETECEPT L001) 

MRITEC (PT L025 

WRITE (LOT 61903) 

ME ITEC TPT. 1004) 

WRITECTPT, 1905) 

RETUIN 

LOOT FRAMAT CAD Ke SI Oke THEY aks OHH Sy FUL Ge FUL! AKG SS 65K 
LY AAABAA*, SK PPPPPPAYLSK, LOD 

LOOL FOPMATES2 Ke SI Gy SHH GOK S THE pK eT UUS POE SUNY OBES // SV 5Ky 
UPAR aK PAA GK PONV GK, IPD HGH, ELOY 

LOO? FORMAT C32 SIT ay THLE HMI aH UU Ry UN TRE TL Oly 
LTABAAABAA TYG, *PPPOPPPY 5K, LL ED 

POOd FORMATO IZ Ke SI! eG THY poh, MHHY, Oy MUU hy Oy CU OR S/he Ky 
te OR PANY OXF PPY LOK LLND 

100% FIPMATE 25K SIG Ke IE Ge THN Ke THM GHG TUT OOK H UU GSK ELLE, 
LRH SAAT GK OT AAT Ore PRE LORS LESD 

LOVS FORMATE QT Ke SILI I IG TH AGTH GS eo SUUJUU! pK SST Ie 
LOAM Ky BAN Ke OP EL OK, LELLELEL DD 

L006 FORMATE STHS LL LIK ef OO CICOt SH, TRAE REAR YOK, AAABAL ST, SKY TYNE OKy 
TENN) 

1007 FCPMATC3FX. LLY 
VOX PMN TOKEN 

LOOR FORMATE ITA, *LL* 610K 
DENNS gL eI NT ZR ye ONE? 

LOOD FORMATE IT Ks LEM LOG IOt pOKe OO! mR eh RAN 2K MABRY AKG HANTS K GT AAN, 
LOK TNE 2K e ING Ly ENN ED 

LOLO FOBMATUSTH SEL! LOK et MON SE POO OK "RAN IK, 2A SKy AAT OK, AAT, 
16K INNES Ky INANE 

LOLL FORMATE AT Ke LUELLELLL 65x, OCOD 
Lat eth NMED 

LOL2 FORMAT (32K "SSSSSS%eGXe TLL TLL Oke 
Leaasease,§: tye 

LOLS FORMATU ALK SSeS Ke She TELL The GG! She! GK 
LOK, f AAT oe BAT some LL? 

VOLS FORMAT (32% PAR LL TRG GGT LOK y PNN YG IEG Nh 2K, 
Leaaaays * 

1015 FORMATUAT KE 
14x, "AF ax 

LOLS FORMATO SLX SS) Tee Itty 
1GX TAA ay of LE | 

LOLT FJQMAT( I2Ke*SSESSS GRIT ILL poke GOGGGG" oS et NN PON yt NN® pay 
TAA GK PAA SK TLELELLEL YD 

VOLA FIRMAT( LK EFEEEEEE Qu yt NNT any! 
USXeMAAROTPA AKON IOI SK, TNA Ky IN’ 
ZPEFEEFE EE Oko MN GK NNT, OR, TET ET OTT!) 

LOUD FORMATCIK TEES LOR GENANE SK GEMNY pe ETM VE GOK PUNE TEGO LES TK 
VGH TRAN GK TOD ehh TIONG Ke INNNE GF IHD PMN ET aE TMUME 2K, eM 
TeEE* tox 
Q*NNN Gy Dey INN Tae OTTO 

LO20 FOOVATU IK PEEFEEET BLINN TS tH, ENS? 

BUR OIREY SKF ON on yt 


LOK eOF eRe (CO Ke MRR OKs (RAH OK® want AAty 


’ 


D+, X_N eX, FREPRRAY -SKeARANALSA pOKy 


Sky MARU Gy ZEN GK, PART GK, TART, 


SGGGS* eSKeINNE EIN Shy 


NO SK EENNT 


"160, 


hal) MC 
UU) 


Pe Me GGG Oe TMNT 2H INE TG EINE 


eS Ke "GG" ORO NNT UES ONENS 


Vi eee ety See CLTET Es 
OKT UM OK UME GY 


soar, 
ake 


He IVOT OR UV GIRO TE Ss 
TOL PNG  2Ey ENN ORG EMM O TKS 


MO we She tNN NNT) THe OTTED 
NNO 2KeI NY ol hy PU ere ee 


ROO By OIG OKs FN y Oke PNY 2K y ING G DEG PNM cu OEE, 


17x, '@2",2K 


Peet Trey 
BES WYNN WRG TT Tee tae aR, 
Nt eGR EMMY OT IMME GES TEE TS LOKy 


tor? FCe@warcix 
leaneysey? 
LINN, 
2d RNN ETH OPTED 

LOZS FORMATE IRA TEEFESEE ET, Gat NN Om 
LO 929 She FOOT ID Ske INN ON 
ZPEREESEEEt OR GIANT oe NMS IRS ITT 

2090 FORMATUIMI S741 

ZOOL FOewAriys) 

2002 FOAMATI ISK SISK saat) 
END 


WoT Rel VV eK oo TELETE 65h e FRI, 
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FORTRAN MOOULF (LEST »NaMS UD 


CARO ND eee CONTENTS meee 


1 SUaEQUT INE DATE( Te IDeIM) 

2 OLSENSECN MOCL2) 

3 Data 23/31 p60e Fle L2VeL520LB2 0213 02449276305 0335 03067 
»* c AY NUMBERS ARE SET £02 LEAP YEAR 
5 IPT=6 

6 00 1b sal, 2 

v {FCL.6T MOtJ)) GO TOL 

8 Imes 

9 to#l 

to LF(J.FV1) GOTO 2 

tt to=1-"0(J-U) 

12 Go Td 2 

1) 1 CONTINUE 

14 2 QETURN 

15 END 
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AUTOFLOW CHART SET = LORAN GROUP PHASE LT A@ARY 


CHAPT TITLE ~ SUBROUTINE DOTIEC AIM SI De TY) 


. 09/28/To INPUT LISTING AUTOFLOm CHART SET - LORAN GROUP PHASE LI ARARY 


gy FORTRAN MODULE (LTST,NAMSQ) 
caso NO eoce CONTENTS eeee 
- t SUBROUTINE OOTOBLAsIMsID,I¥) 
2 OOUALE PRECISION 4e8+C+0 
3 DIMENSION AlOeS dy Sloe LL ICO S40eM12) 
* EQUIVALENCE UBL Deb d eCAULDD oLAEL LOD DELI? 
L OaTa C/* 009000 *+*00 00%, "00 00*,*00 
6 1* 000000 *,* 1 a: 
7 2° SEEPS 22222 * 22s 
8 3°22222222%,* 333333 * 33° 
9 4° 333333" 64at 4% 46" 
to Md 4969555595959 955 . 
1 6" 555555 666666 '4'66 6 
12 T* 666666 *.* 77777777! WwW 
Bb oem *,* BABBRA & 88° 
16 9* 888538 '/ 
1s DATA O/* $999599 "5°99 99% 6! 9999999" 5! 9986999 a9", 
16 B* 999999 "o* gla CE Mee FL” Fee us wee te *e 
7 2° tt nl 
18 K=IM/LO*L 
19 DO Ll {#166 
20 1 ACT sLO=ACT eK) 
2u K=IM-104( K-L)eL 
22 00 2 I#1.6 
2 2 ALT e2V =A E eed 
24 OF 3 1*lee 
25 ACT,60280T 21) 
26 DACT,3=ACTethd 
. K=10/10¢L 
23 00 4 [*le6 
29 @ ACL SSP BCL 
30 K=I0-10°(K-LbeL 
3 DO S [leo 
32 S ACT +S5)=B801 4") 
33 00 6 t#l6 
34 6 ALT + 7)=BCT 8d 
35 Kelyel 
st 00 7 [=leé 
37 TACT 8) =8(T eK) 
38 RETURY 
39 END 
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IN@UT LISTING 
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CONTENTS 


SUAROUTIVE SCOTAC Tay ICN NC ye LARRe Sed 
DIMENSION [aRetyd 

1PT=6 

IC*tin 

Nene 

KO 

TBND#IC eN=1 
TFCES-Twd eA. T.TAND-8) GO TC 900 
NweClANDe 71/78 

OO Lb EwleNw 

LelIC-Lies 

1t=0 

(2=faR2(fwel-t) 

CALL SHFTLOIL E2eb em? 

ICF sfAND 

TFC ICE.GT Bb ICF=8 

OO 2 L=IC,ICcF 

tl=0 

CALL SHETLOILLe L2eoeM) 

KeKelootl 

CONTINUE 

NeN-ICF efC-1 

TBND=N 

{cl 

CONTINUE 

RETURN 

WlITEC 1PT, 1000) 

RETURN 

FORMAT( "CHAR REQUEST EXCEEDS APRAY SIZE*) 
END 
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CONTENTS 


SURROJT INE CLEAR(K YN) 
DIMENSION XKOND 

OO L [#1sN 

x(1)#0.0 

RETURN 

END 
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1 ERwistereiarse ft 


1 ser re 
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CHANGE DDNAME 4 


INCREMENTING SENPUTLE 
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tonne cH . 
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THE SAME PROCEDURE AS 
THOUGH 4N ENO-OF-FILE 
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cette | OL 
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(GET 032801 02 
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Oe eee . 
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@--~~--------------- . 
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oSHETE Os 
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41,068-->1 
| sore on 


eevee eeevee 
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I 


foooo10) 1 02 
a eee 
' SAVE ! 
1 (Loy l2bee® ' 
e-——-—---—_ — —---« 
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' 
coooor2) | 03 
e—-——- = +5 — -—=- . 
1 ua 12,15 ' 
1 
e sk L3,SavEee ft 7 
1 1 
1 ur 2513 i] 
1 ibe 
1 ta 13, SAVE | 
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Bae Soe eg ee eee 
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Appendix B 


VOLTAGE MEASUREMENTS AND OTHER DATA 


This appendix contains a tabulation of the means 
of each of 32 voltage measurements made at each field 
site and simultaneously at NAVOBSY. These means are the 
averages of second-stage edited data taken over 20-min 
periods. At each field site and at NAVOBSY data were 
taken for 1 h. Therefore, the data are conveniently 
broken into three segments, i.e., the first, second, 
and third 20-min periods. 


The mean voltage measurements in volts are 
presented in Table B-l. This table also contains 
"other" data. These data are the year and day number, 


hour of the day, minute, second, block time, time 
difference, distance, and attenuation (see the key 
tabulated on the first page of Table B-1l). Time is 
UTC. The "block time" is the time of the first mean 
voltage measurement in the data segment in milliseconds 
modulo 1 s. The "time difference" is the difference in 
time at the field site relative to NAVOBSY. The 
"distance" is the distance from the measurement antenna 
to the transmitting station in kilometers. The 
"attenuation" is the attenuation in decibels that was 
introduced into the measurement system to achieve the 
desired voltage at the A to D converter (see Volume C). 


The nominal times of the voltage measurements in 
microseconds relative to the start of the pulse are 
given in Table B-2. There are two exceptions. All the 
Measurements made at the field site at Marietta, OH, 
were made 5 us later than the nominal values, and the 
measurements of segment 2 made with the field site 
equipment at NAVOBSY on 1975 day number 188 were made 5 
us earlier. 


All the voltage data of Table B-1 are plotted in 
Figs. B-l through B-12. The sample points where’ the 
Measurements were made were located on the pulse such 
that the odd-numbered samples were near zero crossings 
and the even near peaks. Thus, the odd-numbered 
samples were very nearly in quadrature with the carrier 
and the even in phase (see Fig. 2 in the main body of 
the text). 
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The quality of the data is indicated in Table B-3 

by the signal-to-noise ratios of the edited data. 
These signal-to-noise ratios were computed on a "per 
pulse" basis. The signal-to-noise ratios shown for 
stage 2 edited data in Table B-3 apply to the data 
tabulated in Table B-1 and plotted in Figs. B-l through 
B-12. 
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‘ Table B-2 


Nominal Time of Measurement 
Triggers from Start of Pulse 


Measurement Time 
Number (us) 
1 5 
2 Pas 
3 10 
4 LLY 
5 15 
6 i OST are 
7 20 
8 By Pes 
9 25 
10 7) Me 
11 30 
£2 32.5 
13 35 
14 BF dare 
15 40 
16 42.5 
17 50 
18 OW re, 
19 65 
20 7265 
21 80 
Ze 67.5 
23 95 
' 24 102.5 
25 110 
26 117.5 
27 125 
28 132.5 
29 140 
30 147.5 
31 155 
32 162.5 
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Appendix C 


PROGRAMS FOR ANALYSIS OF DATA 


This appendix presents the computer programs that 
were used in the analysis of the data. These programs 
are written in a computer language known as APL (A 
Programming Language). See Refs. C-1 and C-2 for the 
exact interpretation of the language as executed on the 
Laboratory's IBM 360/91 computer. The program entitled 
TFIT fits a polynomial to the data of Appendix B by the 
least squares criterion. The remaining programs are 
straightforward algebraic manipulation and plotting 
routines. One of these, SLEW, requires elaboration. 


As mentioned in the main text under "Analysis," 
the polynomials fitted to the I and Q data were jointly 
adjusted so that the Q polynomial represented a 
condition of minimum energy. This adjustment is arrived 
at as follows: 


V(t) = I(t) sin wt + Q(t) cos ut, (C-1) 
where V(t) is the voltage function, 
I(t) is the in-phase modulation, 
Q(t) is the in-quadrature modulation, 
w is the angular frequency corresponding to 


100 kHz, and 
t is time. 


The problem is to find an adjustment, At, of the 
sample points in time that results in minimum energy in 


Q(t). 


Introduce At into Eq. (C-l). | 
V(ttdt)el(ttat) sin w (ttdt)+Q(ttdt) cos w (t+dt). (C-2) 
Rewrite Eq. (C-2) as | 
V(ttst) = I(t+At) (sin wt cos wittcos wt sin wit) 

+ Q(tt+dt) (cos wt cos wdt-sin wt sin wit). (C-3) 


Rearranging terms, Eq. (C-3) becomes 


MORE ve ALE 
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V(ttAt) = (I1(ttdAt) cos wht-Q(ttAt)sin wAt) sin wt 
+ (Q(ttAt) cos wAttI(tt+At)sin wht) cos wt. (C-4) 


The energy proportional to that of the cosine 
modulation is observed to be 


totat 
EQ(t + At) = f Q(t + At) cos” wAt dt 
t,tdt 
totat 
+2 i. Q(t + At) I(t + At) cos wht sin wAt dt 
t, tat 
totat 
2 2 
+ I(t + At) sin” wht dt. (C-5) 
t,tdt 
Since the energy in the sine modulation as 


measured is large with respect to the cosine due to the 
locations of the samples on the pulse, it can be 
expected that a solution for At will produce very small 


~ values. By comparison, the interval over which Eq. 
(c-5) is to be integrated is arbitrarily large. 
Therefore, a reasonably accurate solution for At 
corresponding to minimum energy in the cosine 


modulation can be obtained by letting ttAt=t in 
Eq. (C-5), and then taking the partial derivative of 


EQ(t) with respect to At and setting it equal to zero. 
Thus, 


=—255 = 


MA te 


PELL: shad Tins it 


~ tw e " s x 2 
* eh Te Si RS TY a TR ee Be har = 


to oo) 


sin 2 wit f rt) dt - sin 2 wit f Q7(t) dt 


ey ey 
fo 
+ 2 cos 2 wit , O€e) Bde) de = 0 , (€-6) 
fy 
and 
Fo 
2 f Q(t) Le) dt 
t 
pi 
At = x tan t, t, ; nF) 
2 2 
i: Poe) ac - ff TG) lt: 
ty ty 


SLEW executes Eq. (C-7) to find At. 
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APL PROGRAMS 


The program CASES calls the program TFIT to fit 
polynomials to sets of data identified by JMIN through 
JMAX. 


VCASFS(QRIV 

V CASES 
[a J+JIMIN 
pa) +CLFARXiWHICH='Q! 
{3] ICT+((D-1),9)p9 
C4] STARTI*«SDI*«BIGI+MUT+10 
£5] +DO 
C6] CLEAR: 9CT+«((D-1),0)p0 
(7) STARTQ+SDQ+BIGQ+IHUQ+1 0 
[8] DO:70x1J=IMAX 
[9] J+J+1 
E40] FIL 
(11] +SAVQxi.WHICH='0O' 
Eto “SAVEE 


[13]  +DO 
[14] SAVQ:SAVEQ 
15] +b0o 

V 


CASES also calls SAVEQ or SAVEI, depending on 
whether Q or I data are fitted. SAVEQ and SAVEI 
arrange the results of TFIT in suitable arrays. 


VSAVEQLDIV ES 
VY SAVEQ 
bet £4) QCT+QCT,QC 
£2] MUQ*MUQ, 14+MEAN 
C3] BIGQ+BIGQ,~14+MAX 
C4] SDQ*SDQ, 1+SIG 
[5] STARTQ+*STARTQ,TO0 


VSAVFICNIYV 
V SAVET 
baa LOPS CICe LC 
C2) MUI*MUI, 1+MEAN 
e'3) BIGI*BIGI, 14MAX 
C4] SDI+*SDI, 1+SIG 
(5] STARTI*+STARTI,70 
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The program TFIT fits a polynomial by least 


squares to a set of either Q or I data identified by J. 
The degree of the polynomial, D, must be at least 2. 
The first two coefficients of the polynomial (ascending 
order of degree) are set to zero. 


vies Me sce 18 fas 0 set UE gies FA ts Wh ce Dae 
MNIAOANEWNHEH 
Lye tbtsgbsubistiti ti 


v 
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VTFITCOIV 

TET P 

AL+TO+K*LL+NN +0 

L+2 

Net 

MAX«MEAN+SIG+#10 
DO:+*0x1(L>12)4+(N>11) 

K+K+1 

+9x1WHICH='Q! 
IT+0.1x(-70)+(2.5+5%*18), 
42.5415%*18 

IC+IDL; ,JIBITo. *141D 

RES+( (IT. *141D)+.*%,IC)-IPL 3d] 
+MORE 
0:907T+0.1x(-70)+(5x18),35+15*18 
QC+ODL;,J]/POTo. «141d 
RES+((QTo. *141D)+.*%,QC)-9PL 37] 
MORE :MAX*MAX,1A4RESCYV| RES] 
MEAN+MEA2 ,(+/RES)+0RES 
SIG+SIG,((+/(RES-MEANDCK])*2)#((pRES)-D-1))* 
0.5 

+SKIPx 1 K#1 

TO0+TO0+AT+1 

+DO 
SKIP: +(SKIP1AxSIGL K-1xAT>0J]>SIGLK-1xAT<0])+ 
SKIP2xSIG(K-1xAT>0)]<SIG(LK-1xAT<0] 
SKIPI:NeN+LL 

LL+0 

Niv«+2 

TO+TO+AT#122*«L 

+DO 
SKIP2:L«L+NN 

NN+O 

LL«<2 

TO+TO+AT+#-142%N 

+DO 
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The program TRANS! translates the polynomials 
produced by TFIT to rereference time by a time 
increment, TO, which is determined iteratively by TFIT 
to obtain the best fit. 


VTRANS1(O1V 
V TRANS1 
J K+0 
J] DO:70xi1X=ELE 
] K+K+1 

] ITVP(L 3K]+«(-0.1x(VAN/STARTI)(K]) TRP IPVPE :X) 
5) QTVPL sKJ«(-0.1*(VAN/STARTQ)(X]) TRP QTVPL 3X] 
J ITOP( 3K]+(-0.1*(O3S/STARTI)(C¥]) TRP ITOPL;X] 
j QTOPE sKJ]+(-0.1*(O3S/STARTQ)[K]) TRP OTOP[;K] 
J +DO 


The program TRANS2 translates the polynomials 
produced by SLEW to rereference time by the time 
increments, ATV and ATO, which is determined by SLEW to 
achieve minimum energy in Q. 


VIRANS2(01V 
VY TRANS2 

C1] K+0 
£23 DOs+Ox tK=ELE 
[3] KeK+1 
C4] ITVS(;K]<(ATVEX]#10) PRP ILVRE 3%) 
£5 QTVSE ;KJ<CATVLK]210) TRP OTVRE 3X) 
C6] ITOS( sK]«<(ATOLK]#10) TRP TTORL 3K) 
oR) QTOSE ;K1]<(ATOCK]:10) 2FRP 2TORE 3K) 
[8] +DO 
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The program SLEW executes Eq. (C-7) to produce 
polynomials that describe the condition of minimum 
energy in the cosine modulation. 


VSLEWLOIY 
V SLEW 

Ft) Le 1 7 
[2] QCO+«QTOP 
ES] ICO«ITOP 
C4] QCV+<QTVP 
[5] ICV+ITVP 
[6] RLE1+pQC0 
Ere] ROW<14p2C0 
[8] ATH*+THO*+THV+10 
[9] J+0 
(10] DO:+ENDxi.J=ELF 
C11] JS+J+1 
[242] INTI+~L PI IFECOEsd) PP QCOl:7] 
{13} INT2«L PI rcols;J] PP Tcot;sd) 
C14] INT3<L PI QCOL;J] PP ACOL;7] 
C15] THO«TYO,O.5x 30(2*xIV71)+INT3-INT2 
[16] INf1«L PI ICVL;J) PP ACVL[ 32) 
17] INf2+L PI ICV(3J] PP Icvl(;3J] 
(18] JINT3«L PI QCVE;J] PP QCVE[;3J] 
[19] THV+THV,0.5x 30(2xINT1)#INT3-INT2 
{20] +DOd 
[21] BND: QTVR+(QCVx( ROW, ELE )p 20THV)+ICVx( ROW, ELE )p10T"V 
[22] IPVR+(-QCVx(ROW,ELE)p10OTHV)+ICVx(ROW,ELF)p 20THV 


(23) QTOR+(OCOx(ROW,ELE)p 20THO)+ICOx(ROW,ELE)p10THO 

C24] ITOR+«(-9C0Ox(ROW,fLE)p10THO)+TCOx( ROW, ELF )p 20THO 
= [25] AV+THVx10+02 

[26] ATO+THOx10+02 

{(27] Af<ATV-ATo 
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The program ZEROS computes 


crossings in pulses described 
produced by TRANSI1. 


] 

J 

] 

4] 
5) 
(6] 
7] 
8] 
3] 
10] 
aes 
[12] 
(13] 
(14) 
(15) 
(16] 
E27] 
[18] 
(19] 
(20] 
(21] 
[22] 
£23) 
C24) 
(25] 
C26] 
C27) 
[28] 
[29] 
[30] 
[3n} 


VZEROS(OIV 
ZEROS 
T+#2.25+0.25%1321 
PT+T(2)-7(1] 
T+T+10 
RO+RV*+AZTO*AZIV+ZT0+Z2Z7TV +10 
J+0 
DO:+0x1/=36 
J+J+1 
CMO+(((pT),1) 97) 1990707 37] 
SMO«(((o7%),1)o7)1OITOP[ 37] 
CMV<(( (pT), 1)pPTIsOOTVPL 37) 
SMV+«(((o2),1)oLf)sOLTVPL 37) 
S+1002xT 
C+2002xT 
LO+(SxSMO )+CxCH0 
LV+(SxSMV)+CxCMV 
ZCTO*+ZCTV+10 
JJ+0 
DO3:7*SKIP3xi1dd=(p7)-1 
Jd+JJ+1 


the time of the zero 
by the polynomials 


+SKIPx1((xLOCIS])+xLOCdJ+1]) #0 
ZCTO+ZCTO,2.254(JIxDT)+(DIxX LOT II]) + LOC II I-LOCIT +1) 
SKIP :+D03x1((xLV[dd])+xLV[dJ+1]) 70 
ZCTV*+ZCTV,2.254+(ddIxDT) +(DTXLV(i Id) ) FLV ds J-LVi de +1) 


+D0O3 
SKIP3:+SAVFx1(pZCTO)<17 
Z2C0T0«1+2CT0 


SAVE :AZTO*AZTO,ZCTO~-Sx1pZCTO 


Z2T0«ZT0,ZCTO 
AZTV*AcTV,Z2CTV-SxipZecTV 
ZTV*+ZTV,2CTV 

+DO 


The P in lines 8 through 11 was changed to 3 to compute 


crossings in pulses described by polynomials 
produced by TRANS3. 


zero 
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The program SLEWZ produces polynomials that 
describe pulses that have the zero cressing at the end 
of the third cycle aligned with the 30-us sample point. 


VSLEWZ(DI]V 
VY SLEWZ 
ATV+AZTVE 36] 
ATO+AZTOL 36] 
THV*ALVt10202 
THO«+ATO%#10+02 
QCO<QTOP 
ICO+ITOP 
QCV+OTVP 
ICV+«ITVP 
ROW*<14+0QC0 
PLE+1+pQco 
OTVR<(OCVx( ROW, DLE) p20TNYV)+ICVx(ROV,ELZ)p10rTrv 
ITVR«(-OCVxK( ROW FLEE )p1OTHY)+ICVx( ROW HL )p20THV | 
QTOR+(QCOx( ROW ELF )p 20THO)+ICOx( ROW, FLE)p10TN0 
TTOR«(-QCOx( ROW ELE) 1LOTHO)+ICox(2OW*, ELE) p 20710 
AT+ATV-ATO 
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TRANS3 translates the polynomials produced by 
SLEWZ to rereference time by the time increments, ATV 


and ATO. 
~ VTRANS3(OIV 
VY TRANS3 
/ Eas K+0 


[2] DO:+0x1\K=ELE 

[3] Ke +1 

C4] ITV30 3K]<(ATVEQX%IJ210) PRP TTVIL 3K] 
C5] QTV30 sK]<(ATVO[KX3410) TRP QTVRO 34] 
C6] IT030 3Kl]<(ATOCZIJ+10) TRP ITOREL 3K) 
C7] aQ7O3C;X]+(aTolk]+10) PRP QTORE 3X) | 
[8] +DO 
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The program LORAN evaluates the slewed polynomials 
that describe the sine and cosine modulation, 
normalizes the results to unit energy, adjusts the 
attenuation of the received signal by the amount of the 
normalization factor, computes the pulse envelope, and 
plots the results. 


VLORANLCIV 
VY LORAN 

£2) ROT+<O 
£2) T+#0.1*%2.54+2.5%*163 
C3] FQO+((( eT), 1)oT)L$0OTOSL 37) 
C4] FIO+(((pT).1)9?%)1OITOSC 37) 
[5] FOV+(((CpT),1)oT)100TVS[ 37] 
[6] BIV+(((pT),1)pT)1LOITVSL 37) 
[7] NO+#(+/( EQO%2)+FIO%2)*0.5 
[8] NV#(+/(FQV«2)+FIV*2)%0.5 
C9] FQO<FQ0+NO 
[10] FIO«FIO+%O 
C11] MQV+FQVtNV 
C12] FIVeEIViNv 
[13] FO«(( FQO*«2)+2I0*«2)*90.5 
C14] FEV#((FQV*2)+EIV«2)*0.5 
C15] BRELO«+DBOLS]+20x10eNO 
[16] BELV<DBVLJJ]+20x1001V 
ED7a +0] 
£18] OVE: kK+o0 
C19] DO:WAIT<" 
[20] "LORAN-C SIGNAL VOLTAGE PLOTTED AGAINST TIME IN MICROS 
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ECONDS 1975 DAY NO.:', 4 0 WAY[] 
E245 “PERID SkIR: “SST TRhs-1,* STATION: ',STALJ;]," xm 
TR: ',(2 0 #XMTR(J])," DATA SEGMENT: ',SECLZ;] 

[22] ‘ATTENUATION: ‘',(6 3 ¥8FLO),' DB AT WAVOBSY, ',( 
6 3 WBELV),' DB AT FIELD SITF* 

[23] ‘SECONDARY PHASE FACTOR: ',(6 3 ¥SPFLJ7]),' MICPOSFCON 
Ds! 

C24] ¢ 

C25] '. +«+NAVOBSY + ++FPIEFLD STITT! 

C26] WArIT+™m 

[27] +S¥IPxiK=1 

[28] 50 120 FPLOT(FIO,"IV,FO00 AND FOV) VS Tx10 

(29] TWO:K+1 

[30] +DO 

C31] SKIP: 50 120 FPLOT(EO AND EV) VS fx10 


v 
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The program EDGE evaluates the polynomials 
produced by TRANS3 over the leading edge of the pulses, 
normalizes them to unit energy, and plots the leading 
edge of the envelope of the pulses. 


VEDGELNIY 
V EDGE L 
ROT+0 
T*9+176 
EO*+EV<10 
K<0 
DO:72NDx*1K2=3 
Je Lek 
K+X+1 
"3" BVAL fT 
NVE<((LOV«2)+EIV%2)*0.5 
NOE+((2QO*2)+8I0%2)*%0.5 
ZOtEO NOR TIOF[L1+V¥/07 ] 
RV<EV  NVEENVECAUAP¥VVE I 
+D0 
BEND : WAIT CL 
L 1 Git VOLTAGE PLOTTED 
1975 DAY NO.?*, 4 0 ¥DAYEL SI 
e SELES ' eer robes aa STATION: SPATS: 1, 
KMERS * 62 KHIR 1) ,,* DATA SEGISH Ls ALL* 
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' . tAVOBSY + «9PINLD SITE! 
VAIT<M 
Seo 
#O 120 FPLOT((StE0), (StEV), (54(2*S) 420), (S54(2*xS) tEV) , ( 
(-3)4+BO) AND(-S)4#2V) VS 7,7,(104T7),(10+7),(20+7) AND 
20+T7 

q 


EVAL evaluates polynomials identified by A and J 
at times specified by T. 


VEVALCT IV 
4 FVAL T 
TMOO+(((p7),1)07210)1O(2'97T0' AL J] 
EIO«(((p7),1)p7#10)19(e'r70',A)L 37) 
FQV+(((p7),1)0 7819) 1O(2'O7V' ATs] 
ErVe((CoT), Lo Tti0)saoCeti?V' .4yC3c) 
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The program TRP translates any polynomial, P, by a time 
increment, A. 


VTRPLNIYV 


Et] 2+(((1-19P)O((p?),(pP))pAx 1tipP)x(" 1tipP)eo.! 1+1pP)+. xP 


v 


PI integrates any polynomial, P, between limits 
given by the two-element vector, L. 


VPTLOIV 
VY ReL PI P 
bia P+PtipP 
£2) Ret /Px((14L)*19P)-(1tL)*i1pP 


PP multiplies two polynomials, Pl and P2. They 
need not be of the same degree. 


VPPCOIV 
VY P*P1 PP P23;D 
eee P+P10,xP2 
al D#1+oP 
ES] P<(1-1D)OP,((0,0-1)p9) 
C4] P++/({1] P 


These programs require P_ to be arranged in 
ascending order of degree starting with 0 and all terms 
must be accounted for. 
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The program ECM computes the change in ECM, AECM, 
between the field site and NAVOBSY from the polynomials 
produced by TRANS3. 


VECMCCIV 
Vv ECH 


Et) AECIiI+10 

£23 TeG4t+121 
C3] 9 JedMIN 
C4] DO:%Oxr1JEIMAX 
C5] J«J+1 ' 
[6] Eg EVAL # 

E2] NVE*((BROVe2) +R IVe2)*0.5 
C8] NOR#+((8O0%2)+EIO*2)%9.5 f 
[9] NV«tVVELAtPVVE] 


[190] NO+«tNOPLitVvvy0r] 

[ti] "“S" PVAE ‘go 

[12] 3V+VVx((COV+EOV) *2)+(IV+EIV) «2)% 
65.5 

[13] FO0«+NOx(((Q0<EQ0)*2)+(IO+EIO)*2)* 
6.5 

C14] OfODE;J]+1 PP QTO30 3J1xNO 

£15) FrTopl;:7]«1 PD rrTo3l3s7]xKoa 

C16] QPVDE;J]<1 PD OTV3LE 371 *NV 

C17) ITVDC37I]+1 PD ITV3037)]xAV 

Cis] "DD! BVAE 30 

[19] DEV<+((IVxEIV) +QVxEQV)+EV 

C20] D#O*+((IOxEIO)+00xE90) #20 

[21] AECM+AECIH, (BY-F0)%-0.05*xDEV4DEO 

B22) “D0 

] 


PD computes the Dth derivative of the polynomial, 
P. P must be arranged in ascending order of degree 
starting with zero and all terms must be accounted for. 


VPDLTIV 

V RD PD P3kK3b 
Lig K+0 
EZ] Lep? 
C3] DO:+0x1.K=DLE 
C4] ReP414Px 1419? 
bo K+K+1 
C6] +DO 
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The program DBEL computes the adjusted attenuation 
that is plotted in Fig. F-16 of Appendix F. 


VDBEL(OIV 

Vv DBFL 
£2] M+10 
2] LOSS*ATTENL ;2]+1 
[3] NDB+LOSS(112] 
C4] K+0 
Go J+13 
C6] DO:+ENDx\K=2 
EPI K+K+1 
[8] N+20x10@DIST(JJ#DIST([1] 
C9] NN+LOSS(11-N 
(10] M+M,NNN+NN-LOSS([J] 
[11] WDB+NDB,LOSS(CJ+0,18]+NNN 
[12] J+22 
[13] +D0 
[14] FND:N++/NtpN+LOSS(30+16],ATTFNE(112);1],ATTFNE (304+16)3;1] 
[15] NWN++/NNtoNN+ATTEN( (12+19)31] 
[16] WNN++/NNNtoNNN*ATTENE (21+19)31] 
[17] NDB+NDB,(3pN),(3p9NN+M[1]),3p9NNN+Mé[2] 
[18] DIST+DIST[130],(3pCARB) ,(3pDANA),3pNANT 


= 221 = 


+ AEA ee 


oe 


Pee ae 


en eae encanta NBS NG TLS ESE ITA ENO U oho , 


REFERENCES 
C-1 Anon., "APL Language," IBM Corporation Publication 
GC26-3847-0, March 1975. 


C-2 Anon., “Graphs and Histograms in APL/360-370," IBM 
Corporation Publication S$B21-0415-0, 1972. 
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Appendix D 


POLYNOMIAL COEFFICIENTS 


This appendix records the values of the 
coefficients of the polynomials that were used in the 
analysis of data. Tables D-1 through D-12 tabulate 
these coefficients. The coefficients are presented in 
ascending order of degree from the top to the bottom of 
each column of eleven coefficients. The groups of three 
columns’ under the name of each field site are 
coefficients for data segments 1, 2, and 3 from left to 
right. 


The tables contain the following labels: 


QTOP and QTVP are the labels for 10th degree 
polynomials fitted to Q data taken at 
NAVOBSY and data taken simultaneously at 
the field site, i.e., the van. P denotes 
polynomials of more permanent interest; 
i.e., they represent the data as taken. 


ITOP and ITVP are defined the same, except the data 
fitted are I data. 


QTOS, QTVS 
ITOS, and ITVS are the labels for polynomials slewed to 
achieve minimum energy in Q. 


QTO3, QIV3, 

IT03, and ITV3 are the labels for polynomials slewed so 
that the third-cycle zero crossings 
align at 30 us. 
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Polynomial Coefficients for Q Data Taken at NAVOBSY 
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Table D-8 
Polynomial Coefficients for | Data Taken at Field Sites 
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Appendix E 


CARRIER MODULATIONS AND PULSE ENVELOPES 


This appendix presents graphical representations 
of the measured pulses. The plots are constructed from 
values of voltage evaluated from the polynomials of 
Tables D-5 through D-8. Figures E-l through E-36 
present the sine and cosine modulations for all the 
sites and all data segments. Figures E-37 through E-72 
present the corresponding pulse envelopes. The voltages 
plotted on these figures are those associated with 
pulses that have the same energy. Tables E-1 and E-2 
present the times of the zero crossings in the pulses 
relative to the reference time of the voltage samples. 
The vertical groups of three rows are for the three 
data segments. The labels Cl, C2 ... CO on the left 
denote the location of the instrumented van during data 
recordings (see Table E-3). 
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Table E-3 


Key to Labels That Denote 
the Location of the Instrumented Van 


co 


Cl 


C2 


C3 


C4 


DO 


D1 


D2 


D3 


NO 


Nl 


N2 


N3 


NAVOBSY 
Wilmington, NC 
Emporia, VA 
Towanda, PA 
Dexter, NY 
NAVOBSY 
Danville, IN 
Marietta, OH 
Georgetown, DE 
NAVOBSY 

Toms River, NJ 
Grottoes, VA 


Bluefield, WV 
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Appendix F 


RESULTS OF PULSE ALIGNMENT 


This appendix assembles the illustrative material 
associated with pulse alignment at the llth sample 
point on the pulse, i.e., nominally 30 js from the 
start of the pulses. Figures F-1l through F-12 present 
the leading edges of pulses so aligned. The voltages 
plotted in these figures result from evaluations of the 
polynomials of Tables D-9 through D-12 and normalized 
to a peak value of l. 


From these values of voltage, ECM was calculated 
as explained in the main body of text. The result of 
this calculation is plotted against SPM in Fig. F-13. 
SPM and ECM are plotted against SPM in Fig. F-13. SPM 
and ECM are plotted against the distance from NAVOBSY 
to the field site in Figs. F-14 and F-15. Figure F-16 
presents the observed attenuation of the signal at each 
site relative to the field strength of Carolina Beach 
observed at Wilmington, NC. Tables F-l and F-2 tabulate 
the times of zero crossings relative to the time of the 
llth sample point. The vertical groups of three rows 
are for the three data segments. Table F-3 records the 
values of the distance from the transmitting station to 
the site, ECM and SPM plotted in Figs. F-13, F-14, and 
F-15. Table F-4 records the values of attenuation 
plotted in Fig. F-16. The labels Cl, C2 ... NO in these 
tables refer to the location of the van during data 
recordings (see Table E-3 of Appendix E). 
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Fig. F-13 The Relationship Between SPM and ECM 
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Table F-3 


Values of Distance, SPM, and ECM Plotted 
in Figs. F-13, F-14, and F-15 


1975 NAVOBSY to 
Day Field Transmitting Field Site SPM ECM 
No. Site Station (m) (us) (us) 
177 Wilmington, NC Carolina Beach “520672 ~3.600 “1.653 
“520672 ~3.600 “1.656 
“520672 “3.600 “2 #793 
175 Emporia, VA Carolina Beach “259900 “1,490 “0.772 
~250900 “1.491 “0-829 
“250900 “1,489 “O,, 832 
183 Towanda, PA Carolina Beach 323487 AetRS “Os223 
323487 2.736 “0.049 
323487 2.720 “9.032 
181 Dexter, NY Carolina Beach 569702 3.894 “0.105 
569702 3.901 “0.355 
569702 3.904 “0.369 
163 Danville, IN Dana “827137 50 562 5 SOL 
"8ig1a7 pe ee "5...867 
817137 “$558 “S32 306 
161 Marietta, OH Dana “376178 ee ee ewe ae 
“376178 “3u065 “eC eteo 
“376178 ip ra Be ge Phe) 
156 Georgetown, DE Dana 150434 0.677 a. 759 
150434 0.680 aete 
150434 0.679 "se 223 
149 Toms River, NJ Nantucket 7273230 72 «O92 2,053 
7273230 “2.974 “1.014 
“273290 “2.084 bl 08 
~ 
170 Grottoes, VA Nantucket 170013 1.891 “0.689 
170013 1.896 “0.985 
170013 1.914 Shao 
168 Bluefield, WV Nantucket 405849 3.019 bowre 
405849 3.023 1.897 
405849 3.031 1.879 
154 NAVOBSY Carolina Beach 0 0.026 “0.010 
0 0.037 0.149 
0 0.044 0.153 
188 NAVOBSY Carolina Beach 0 0.046 0.053 
0 0.048 0.048 
‘ 
2 
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Table F-4 


Adjusted Values of Attenuation Plotted 
in Fig. F-16 


Distance 


Station Adjusted Measured 
to Site Attenuation Attenuation 
Site (m) (dB) (dB) 

Cl 23572 0.9000 43.762 
“0.005 43.757 

“0.075 43.687 

c2 293344 21 584 22.178 
“21.964 21.798 

“21.841 pak Bese ak 

c3 867731 “40.541 3.222 
“40.390 i 

40.340 3.422 

C4 1113946 "47.354 3592 
"47.442 ~3.680 

“47.431 3.669 

D1 85946 ~Pigest 32.526 
“42.255 32.547 

11.199 32.563 

D2 526905 ~ 30.403 13.359 
~ 30.384 13.378 

30.391 '3.372 

D3 1053517 “44,861 “1.099 
“44,790 “0.938 

“44.940 ee 

Nl 384346 724.246 19.516 
23.525 20.187 

23.467 20.295 

N2 827590 “39.701 4.061 
39.992 4.170 

39.695 4.067 

N3 1063425 742.921 0.841 
~ 42.996 0.766 

42.867 0.895 

co 544244 “29.140 14.622 
DO 903083 “42.485 Lz? 
NO 657576 “30.380 13.382 
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Appendix G 


RESULTS OF FREQUENCY FILTERING 


This appendix assembles the histograms and_ the 
plotted Fourier coefficients that were used in an 
attempt to understand the data taken at Georgetown, DE, 
and simultaneously at NAVOBSY. 


Figures G-l through G-8 are derived from data 
taken at Danville, IN, and NAVOBSY during the first 
20-min recording. Figures G-l and G-2 are histograms of 
first-stage edited data taken at NAVOBSY and Danville, 
respectively. Figures G-3 and G-4 present the Fourier 
coefficients derived from Q and I data, respectively, 
taken at NAVOBSY simultaneously with recording at 
Danville. The Danville coefficients are shown in Figs. 
G-5 and G-6. Figures G-7 and G-8 are the histograms of 
frequency filtered data taken at NAVOBSY and Danville, 
respectively. 


Figures G-9 through G-12 present histograms and 
Fourier coefficients derived from data taken at NAVOBSY 
during recordings at Georgetown, DE, during the first 
20-min data recording, i.e., data segment 1. Figure G-9 
is the histogram of first-stage edited data. Figures 
G-10 and G-1ll present the Fourier coefficients for Q 
and I data, respectively. Figure G-12 is the histogram 
of the frequency filtered data. 


Figures G-13 through G-16 are the same as Figs. 
G-9 through G-12 but for data segment 2, and Figures 
G-17 through G-20 are for data segment 3. 


Figures G-21 through G-32 are the same as Figs. 


G-9 through G-20 except that the data were recorded at 
Danville, IN. 
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Fig. G-1 Histogram of First-Stage Edited Data Taken at NAVOBSY During Recording 
at Danville, IN; Data Segment 1 
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Fig. G-2  Histograrn of First-Stage Edited Data Taken at Danville, IN; Data Segment 1 
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Fig. G-3 Power Spectral Density Showing Interference in Q Data Taken at 


Data Segment 1 


NAVOBSY During Recording at Danville, IN 
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Fig. G-4 Power Spectral Density Showing Interference in | Data Taken at 


Data Segment 1 


NAVOBSY During Recording at Danville, IN 


345 


% 


. 


a hie oe 


Sai Sale dee 


EE A BS We el, TA a am 


‘3 


oa” 


BLANK-NOT FI 


7 
& 


RECEDING PAG 


‘ 
| 


000070 = ONS ¥3MOd 9G Wil 


s 
8 
i 
9 
s 
* 
€ 
2 
1 


‘Aen 


SE nee ee 


0000°O = ¥3K0d DV SAV 


7600°0O = B3M0d DV Widd 


€0-GI9T9*O= IN WA KvW 


%0-GOE9T*O= BNIVA NIW 
. sz 

2 Ssz 

© sz 

© sz 

ezSz 

els? 

* 362 

= 69% 

. Boz 

. tz 

s 992 

= $97 

e 9%2 

. £2 

. zz 

. 192 

s 092 

= 6€2 

*8ez 

al ez 

. 9e2 

* Sez 

e 92 

= £2 

. zez 

2 1e2 

° o€z 

bd 622 

. B27 

° i2z 

. 922 

$22 
422 
F272 
z22 
ad Tez 

. O22 

© 612 

. B12 

* £12 

2 912 

asiz 

. wl2 

. £12 

e C3z 

© 112 

* O12 

al 602 

s 302 
bd Lez 

2 902 

s soz 
* 902 

* €02 

bd 702 
bd 102 

* 002 

° 661 

@ 861 

. 161 

2 961 
e561 

2 561 

° fol 
e 261 

. T6T 

. 061 


347 


Data Segment 1 


‘ 


Fig. G-5 Power Spectral Density Showing Interference in Q Data Taken at 
Danville, IN 
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Fig. G-6 Power Spectral Density Showing Interference in | Data Taken at 
Danville, IN 
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Fig. G-7_ Histogram of First-Stage Edited Data Taken at NAVOBSY During 
Recording at Danville, IN, After Frequency Filtering; Data Segment 1 
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Segment 3 
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After Frequency Filtering; Data Segment 1 
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Fig. G-28 Histogram of First-Stage Edited Data Taken at Georgetown, DE, 
After Frequency Filtering; Data Segment 2 
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Fig. G-32 Histogram of First-Stage Edited Data Taken at Georgetown, DE, 
After Frequency Filtering; Data Segment 3 
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